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Abstract 

This report is concerned with describing the NASA Pioneer-Venus 
Reverse Playbach Telemetry Program. This program is a software package 
developed to decode telemetry data received from the Pioneer-Venus Multi- 
probe Mission. The program processes recorded data in an off-line mode 
of operation. The program reads a digital tape containing receiver sof^- 
decisions recorded during the off-line processing, in either the forward or 
reverse direction, of the Pioneer-Venus pre-detection recordings. The 24- 
bit frame synchronisation word of each telemetry frame is found. Next, the 
data are decoded on a frame-by-frame basis using either a sequential decoder 
or a quick- look decoding procedure. Decoded data ore available in printout 
form or may be written unto another digital tape in thi ^roper time sequence. 

This report is organised into three basic categories. First, a detail- 
ed description of installation dependent software and program-user inter- 
action is provided. Second, a discussion of the frame synchronisation 
algorithm and decoding procedures is presented. Finally, computer program 
details ore included in the Appendices. 
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1.0 Ia,troductioa 

The NASA Pioneer-Venus Reverse Playback Telemetry Program is a software 
package which has been developed to decode telemetry data received from the 
Pioneer-Ven\is Multi-probe Mission. The program has been designed to function 
on the computer facilities of the NASA Ames Research Center. It will process 
recorded data in an off-line mode of operation. Specifically, the program 
will read a digital tape containing; 1. ) receiver soft-decisions from the Deep 
Space Network's Symbol Synchroniser which will be recorded during the off-line 
processing (in both the forward and reverse direction) of the Pioneer-Venus 
pre-detection recordings; 2.) time tags associated with the data, and 3.) an 
identification file. The 2it-bit frame sync word of each telemetry frame is 
found, and the direction of the symbols associated with each frame sync word 
is then reversed, if necessary. Next, the data are decoded on a frams-by- frame 
basis using either a sequential decoder or a quick- look decoding procedure. 
Finally, the results are available in printout form or may be written out onto 
another digital tape in the proper time sequence. 

The primary objective of the Reverse Playback Telemetry Program is to 
minimize information loss due to ^aps in the received data. The telemetry 
from each probe will be sequentially -decoded by the Telemetry Processor Assembly 
during the real-time data pi^cessing. These data will contain gaps that appear 
during the time that there is a need to synchronize or resynchronize any of the 
units comprising the telemetry system for that particular probe. These data 
gaps will occur for each probe during the mission at 1. ) the initial acquisition 
of received data, 2.) communications blackout when the spacecraft enters the 
Venus atmosphere, 3.) the time the probes change data rate, and U. ) times 
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when the receivers lose lock due to probe motion caused by Venus atmospheric 
wind ffusts. 

The data which are lost due solely to the synchronisation process will be 
recovered by off-line data processing of the pre-detection recordings made 
from the open-loop receivers. Figure 1 illustrates the Ground Station con- 
figuration. Processing to recover such data will consist of 1. ) playing the 
analog tape of receiver outputs in the reverse direction, 2. ) up-converting 
the signal from the tape recorder to S-band and feeding it into one of the 
closed-loop receivers in the real-time system, 3. ) processing the signal 
throxigh the receiver, subcarrier demodulator, and symbol synchronizer, and 
U. ) recording the soft-decisions out of the symbol synchronizer on digital 
tape. This digital tape is used as input to the Reverse Playback Telemetry 
Program. Synchronization is established for each frame and the data is de- 
coded using the software sequential decoder employing the Fano Algorithm. 
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2.0 Program Usage; Initialization of Program at NASA Aaes 

This section explains the installation dependent sofware developed 
for the NASA Ames Research Center's ] 3 M 360/67 computer system, as well as 
the procedures for utilizing the Reverse Playback Telemetry Program. Two 
procedure definitions (PROCDEPS) have been prepared to perform the tasks of 
tape mounting, storage allocation, and logical unit assignments. This 
obviates the need for the program user to type in the sequsnce of command 
strings which would noimally be required. 

Before exercising the program, the user must first submit both the input 
data tape and the output tape (if one is desired) to the Computer Operator. 
The tapes and a TSS Job Card are submitted to the attendent in the Compu- 
tation Center's l/O Room. Each magnetic tape must be labeled according to 
the recommended foimat shown below. 


User 

"T- 

1 

1 

1 

Volume • 

Mail 

Name 

1 

1 

-JL- 

1 

Name 1 

L_ 

Stop 


•Taue Label 


User Name ; This block contains the name of the user authorized 

to use the computer account along with any .tape reference 

infoimation provided by the user. 

Volume Name; JPL (for input tapes at 16OO BPI). 

REVIDR (for output tapes at 8OO BPI). 

Mail Stop ; This information tells the l/O Personnel where to 

return the tapes when program execution is finished. 


Once the tapes have been submitted, the user must log on the NASA Anes 
TSS system. The LOGON procedure for a terminal (ASCII-type) will be given, 
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However, alternate LOGON procedures can be found in NASA Ames Research Center - 
A Guide to TSS/36 o « Connecting the terminal with the system is accomplished 
by initializing the Ames Dataswitch unit associated with the terminal. The 
system will respond with a hyphen; typing TSS and then pressing the Return 
Key connects the terminal with the system allowing the LOGON command to be 
entered - 


. . . NASA AMES SWITCH A . . . 
-TSS 

/LOGON /userid, jterminalid 
TSS 9.0 PTFU9X 
PI=.8T9 

ENTER PASSWORD 


(system) 

(user types TSS, hit Return) 
(user types in command, 
hit Control S to send 
this and all future 
responses. ) 


where, userid - is the user identification that was assigned to the computer 
account . 

terminalid - the terminal identification number of the conversational 
terminal being used for the task. 

- enter user password after being prompted. 

Having logged on to the TSS/ 36 O system t;ie Reverse Playback Telemetry 
Program can now be Initialized by entering the appropriate PROGDEF command 
(RUNl or RUN2). The first PROCDEP, RUNl, is invoked when the user desires 
only printed output. The second PROCDEP, RUN2, is invoked in the event the 
user has an additional requirement for tape output of decoded data. List- 
ings of the PROCDEP files follow along with an explanation of their operation. 
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RUIT 1 

This procedure definition consists of the following system coicmands - 



0000 

PROCDEF RUMl 


0100 

PARAM $TAPE 


0200 

DDEF FTO3F0Ol,VS,LPRHIT.OUT 

V 

0300 

DEFAULT SYSINX-E 


oUoo 

MTMSG 


0500 

please OBTAIM 9-TRACK, 1600-BPI TAPE $TAPE WITH NC RING 

\ 

0600 

DDEF FT0UF001,PS, INPUT, UKIT«(TA, 9D3) ,DISP«0LD,- 

0700 

LABEL«( 1 ,HL ) , VOLUME* ( , $TAPE ) ,PR0TECT»Y , - 


0800 

DCB- ( DEVD»IA ,DEN- 3 ,RECFM«F , LRECL*10U2 ) 

] 

0900 

DDEF FT06F001,V1,TEMP,RET»T,DISP«NEW,- 


1000 

DCB* ( RECFM-F , LREC^»- i.0U2 ) 

T 

1100 

DEFAULT SYSIHX-G 

i 

f 

1200 

DRIVER$$ 

t 

1300 

PRINT LPRINT. OUT, PRTSP-EDIT, ERASE*! 


lUOO 

RELEASE FTOUFOOl 


1500 

RELEASE FT06F001 

4 

l6oo 

DELETE INPUT 

1700 

ERASE TEMP 


I To begin execution of the program with the printout-only option the user j 

I issues the command; RUNl J?L where, RURl is the name of the above procedure 

command file and JPL is the volume or name of the input ’tape containing the 
data to be processed. It is possible to specify tape names other than JPL 

, i 

with this procedure definition since the tape name is an input parameter 

specified by the user. | 

i 

The basic operation of the procedure command file is to alert the computer i 

i 

■ I 

I operator as to which tape is to be mounted, assign logical Fortran 1/0 devices j 

I 

' for program operation, and initi e program execution. The input tape defin- ]\ 

ition command defines the DDNAME "FTOUPOOI" and associates it with the physical 
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sequential (PS) data set "IITPUT". This data set has been previously written 
(old), the tape is unlabeled, and is to be read from a 9-track tape drive 
having a density of l600 BPI. The volume name supplied by the user is sub- 
stituted in for the dummy parameter and the particular tape is to be 

write protected when mounted. The Data Control Block indicates that the 
physical device is a tape drive (DEVD»TA), the density is l600 BPI (DEN»3), 
the record format is of fixed-length (RECFM»F), and the record length is 10U2 
bytes (LRECL“10h2) . Lines 0900 and 1000 are the deta definition for the 
mass storage. The DDHAME 'FTO6POOI'' is associated with the virtxial indexed (Vl) 
data set named "TEMP". This assigns Poirtran logical unit 6 to the temporary 
disk storage file. The file is specified as being temporary (RET»T) and that 
it is to be created (DISP»HEW). The character of the data set is described 
as being of fixed record forme.c with a record length of I0U2 bytes. Program 
execution is initiated by issuing the command DRIVER$$. The user now will be 
quizzed by the program as to which options are desired. When program execution 
is terminated the print file will be released to the line printer if output 
to the line printer was specified. Finally, the device assignments are re- 
leased and the temporary mass storage file is erased from the user's library. 

Rmr 2 

This procedure definition consists of the following system commands- 

0000 PROCDEP RUN2 

0100 PARAM $TAPEI,$TAPE0 

0200 DDEF FT03F001,VS,LPRIUT.0UT 

0300 DEFAULT SYSDrX=E 

OUOO MTMSG 

0500 PLEASE OBTAIN 9-TRACK, I6OO-BPI TA’’E $TAPEI WITH NO RI1?G 

0600 MTMSG 

0700 PLEASE OBTADI 9-TRACK, 300-BPI TAPE $TA?E0 WITH RTTJTG IN 
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* 0800 DDEF FTOUFOOl, PS, INPUT, UNIT«(T!A,9D3), DISP»0LD,- 

0900 LABEL- ( 1 ,HL ) , VOLUME- ( , $TAPEI ) ,PROTECT»Y , - 

I 1000 DCB-(DEVD«TA,DEN»3,RECFM-P,LRECL-10U2) 

1100 DDEF FT05F001,PS,DATOUTE1,UNIT-(TA,9D2),DISP-NEW,- 

, 1200 LABEL-(1,NL),V0LUME»(,$TAPE0),PR0TECT-N,- 

1300 DOB- ( DEVTD-TA , DEN-2 ,RECFM-F , IRECL-1216 ) 

lUOO DDEF FT05F002,PS,DAT0UTE2,UNIT-(TA,9D2),DISP»NIW,- 

1500 LABEL- ( 2 , NL ) , VOLUME- ( , $TAPE0 ) , PROTECT-N , - 

1600 DOB- (DEVD-TA, DEN-2 ,RECFM-F,LRECL-1216) 

• 1700 DDEF FT06F001,VI,TEMP,RET-T,DISP-NEW,- 

' 1800 DCB-(RECFM«F,LRECL-10U2) 

, 1900 DEFAULT SYSm-G 

I 2000 DRI7ER$$ 

2100 PRINT LPRIira.QUT.PRTSP-roiT.MASE-Y 

2200 RELEASE FTOUfOOI 

J 

2300 RELEASE FT05F001 

' 2400 RELEASE FT05F002 

' 2500 RELEASE FTO6FOOI 

■ 2600 delete input 

2700 ERASE TEMP 

To begin execirtion of the program with the additional capability of writing 
I 

an output digital tape, the user issues the command: RUN2 J?L,RE/IDR where, RUN2 
^ is the name of the above procedure command file, J^L is the name of the input 

tape, and REVIDR is the name cf the output tape to be written. Again, it is 
possible to specify different tape names as desired since the tape names are 
input parameters issued by the user. 

The operation of this procedure command is identical to that of RUNl with 
tha addition of an output tape definition. The DDNAME"FT05F001" is associated 
with the data set "DATOUTEl”. This data set is to be created (NEW), the tape is 
unlabeled, and to be written on a 9-track tape drive ha’rt.ng a density of 800 BPI. 


llie record format is fixed with a record length of 1216 bytes. A double end- 
of-file is written after the last data record on the output tape. 
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A flowcha-rt which illustrates the basic operation of the program is 
provided in Fig. 2. The program begins by setting the logical unit numbers 
for each of the inputs and outputs used by the program. These assignments 
are given in Table 1. These logical unit numbers must be assigned to 
specific devices by the user before the program is executed. These assign- 
ments were discussed in Sec. 2.0. 

The next step in the program is to interact with the user to input the 
various options and values necessary to guide program execution. The options 
the user must specify are asked for by the program in a conversational Kode 
of operation. Specifically, the following is asked of the user; 

A. ) IS DATA HECOBDED ON TAPE IN FOEWAED OR REVERSE DIRECTION (FW^ OR 

REV)? 

The user responds with FWD if data on the tape are recorded in the 
forward direction, or REV if data on the tape are recorded in the 
reverse direction. 

3.) DO YOU WANT TAPE OUTPUT OR PRINT OUTPUT (TAPS OR PRUT )? 

If the user responds with TAPE, program output is written to 
magnetic tape, whereas, if the user responds with PRHT, program 
output will be printed on a hard-copy device. 

C.) DO YOU WANT PRINTOUT TO THE TERMINAL OR TO THE LINE PRINTER (TERM 

OR LPRNT)? 

The user may choose to have the program output printed on the 
terminal for short sessions or sent to the line printer when 
large anounts of data are to be processed. 
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^iHICH PHASE SEFEBENCE DO YOU WAMT - DATA OR DATA-BAR (D OR DB)? 

The user must Indicate the desired phase reference to be used in 
mapping the soft -decisions (see Pig, 10), The user responds with 
D if the DATA phase reference is desired or with DB for the DATA- 
BAR which is a l 80 ° phase reversed assignment. Pano decoder 
operation is sensitive to the phase reference used, hence the 
correct phase assignment is the one which results in a frame being 
successfully decoded when reliable input data ore used.i.e,, a 
frame with high signal-to- noise ratio. 

IS THE DATA TO BE DECODED OR IS RAW SYMBOL OUTPUT DESIRED (DEC OR RAW)? 
The user must Indicate whether data are to be decoded or only the 
raw data symbols (undecoded) are to be printed out. The raw data 
symbols are printed out in octal and binary fora. This question 
is not asked if the user has previously specified tape output since 
raw symbols cannot bo outputted to tape. 

DO YOU WANT THE PANO DECODER USED OR THE QUICK-LOOK OUTPUTS (F.AN‘0 
OR 5L)? 

If the data ore to be decoded, the user must specify whether the 
Pano sequential decoding algorithm or the "Quick-Look” algorit.hm 
is to be used. This question is not asked if the user has previously 
specified raw data symbol output. 

DO YOU WANT ALL THE FRAMES PRINTED OUT (Y OR N)? 

The user must specify whether all fromes are to be processed, or 
if only selected frames are to be processed. Responding with Y 
will result in all data frames being processed according to the 
options thus far chosen. If the user wishes to select the frames 
to be processed a response of N s.hould be issued. 
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:o YOU WAS? TO sracirr ?m frames tc be procxssxc accctoeio to t»<e 
:ao an rum sukber (time or numb)? 

If Zii.9 ustr fptcifics that only Mltetod fraot ar« to bo proooosod, 
it suit bo indlcotod wnotbor tbo frMios to bo toloetod will bo 


rpooiflod by frwo ovabor or tiao tof. If tbo viaor roopondo vltb 
TIME* tbo progroB vlll roply vltb tbo foUowlnf* 




SRCZTY THE TIME TAO Of EACl ERAMB WHICH IS TO BE PRIIRID OUT. 

ISDICATE YOU ARE FUnSRED BY EHTSRnO A *1.0 

HOURS- 

Tbo ubor tbon ontoro tbo hour froB O.C to 23.0 of tbo tiao tog 
for tbo doolrod fr«o. Tbo nuBbor ontorod nuat bo is floating* 
point fozBot. Tbo prograi tbon print i out: 

KUfms- 

Tbo uaor ontoii tbo nlnutoi frcB 0.0 to 59.0 of tbo tlao tag 
for too doiirod froBO. Tbo nusbor ontorod auat bo in floating* 
point fozvat. Tbo progra vill tbon print out: 

SECONDS- 

Tbo uaor ontors tbo aoconds fros 0.0 to 5^.C in floating- 
point fomat. Aftor tbo uaor naa ontorod tbo boura. ainutoa. and 
sococda of tbo tiao tag for a doairod frno, tbo prcgm will con- 
tinually roquoat tbia typo of infcmation for aoro fraoa that aro 
to bo aoloctod by tiao tag. Tbo uaor indicatoa that ro aoro fraoa 
aro to bo procoaaod by onttring a -1.2 aa tbo tiao input. 

If tbo uaor roaponda to quoation H with NUMB, tbo prograa will roply 
wltb tbo following. 

SPSCI7Y BY NUMBER THE FRAMES WHICH ARE TO BE PROCESSED. THE LAST 
Disir cr EACH TJMBER MUST END 'JP 3 CCLLT'W f. 3DICATE YOU ARE DCHE 


BY HTTTrn; A RETURN. 
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I. ) DO YOU WAUT SAW SYMBOLS TO BE SYWCHROHIZED (102U) IN A DATA FRAME OR 

THOSE (102J) CORRESPOHDINQ TO A TAPE RECORD (SYNC OR TREC)? 

If the user req.ue8ts SYNC the row symbols will be synchronized to agree 
with a data frame, otherwise the raw symbols will be printed as they 
appear in a tape record. 

The user must now enter the numbers of the frames which are to 
be processed. There is to be one number per line with each number 
ending in column 5« Frame number input is terminated by either 
entering a zero or hitting a return. 

After having selected the desired options, the user is asked to specify 
input values for several parameters. The parameters axe; 

A. ) The DSS number. ENTER DSS NUMBER (Ik OR U3)- 

The user must respond with either 14 or 43. 

B. ) Probe identification. ENTER PROBE ID (SP1,SP2,SP3, OR LP)- 

The user should respond with SPl, SP2, SP3, or L? corresponding 
to small probe ^ 1, small probe / 2, small probe 4 3, or the large 
probe, respectively. 

C. ) Tape sequence number. ENTER TAPE SEQUENCE NUMBER- 

If tape output was requested, the user must specify a two digit 
tape sequence number. 

D. ) Year of data recording. ENTER YEAR DATA WAS RECORDED (E.G. 78)- 

If tape output was requested, the user is asked to input the last 
two digits of the year in which data was recorded. 
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E,) Decoder computations limit. INPUT MAXIMUM NUMBER OF COMPUTATIONS 
PER FRAME. THE LAST DIGIT OF THE NUMBER MUST END UP IN COLUMN 8. 

0 INDICATES DEFAULT ■ 100,000 COMPUTATIONS. 

If the Fano sequential decoder 'tf«4S requested, the user must specify 
the maximum number of computations allowed per frame (see Sec. 5 . 3 ) 
in integer format with the last digit in column 8. When a zero or 
blanlt is entered, the default value of 100,000 computations per 
frame is set. 

After all the options and parameters have been specified, the contents 
of the input tape are read onto the direct-access device assigned to logical 
unit 6. The space allotted must be sufficient to accommodate a maximum of 
5000 records, each record being 521-18 bit words. 

After the data are read onto the direct access device, the data are 
processed and outputted according to the user's specifications, l^on com- 
pleting the assigned task, the program asks ARE YOU DONE (Y OR N)? If the 
user is finished a Y should be entered to tenainate execution of 
th2 program. Should the user wish to process the data for a different set 
of options, responding with an N is required. Now the program will request 
a new set of options by proceeding as described before. It should be noted 
that if N is the response, the user will not be given the option of choosing 
a different tape direction. This is due to the fact that the program accesses 
the data directly from the direct-access file for the subsequent program 
passes. 

Pig. 3 shows a sample session using the Reverse Flaybeck Telemetry 
Program. Here, print output had been sent to the terminal instead of t.he 
line printer. Note that if the user mistypes a response, the program asks 
the some question again and waits for a correcc response. 
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The encoder aboard the spacecraft Is a convolutional encoder as 
illustrated In Fig.U. The convolutional code employed is non-systematlc, 
i.e., the input data do not appear explicitly in the encoded output sequence. 
Systematic convolutional codes, on the other hand, are characterized by the 
fact that the first stage of the shift register is the only connection to 
one of the channel output ports. The code rate, R, is one-half. That is, 
for every input information bit there are tvo channel symbols at the output 
of the encoder. The encoder has a constraint length of 32 bits. A telemetry 
frame consists of 512 infoimation bits (1024 channel symbols). The last 24 
bits is a Frame Synchronization Word. This is a fixed binary sequence which 
appears at the end of each frame. It has the foim 1111100011000 
1010100100 1 . 

The encoder in Fig. 4 outputs two channel symbols denoted by P and Q, 
where ^ is the complement of Q. The value of each symbol is based on the 
values of 32 selected data bits previously fed into the shift register. 
Modulo-two adders are connected to the various shift register cells as shovm. 
Each symbol, ? and Q, is a logical 'one' if there is an odd number of 'ones' 
in the selected data bits and a logical 'zero' if there are an even number. 
The encoding cycle begins at the end of the last bit of each frame syn- 
chronization word. At this time the shift register and the flip-flop used 
to generate the code are reset to a logical 'zero'. If the inverter on the 
Q channel output port is removed the resulting structure of the set of all 
outpuz codewords y can be viewed in terns of a code tree as shown in Fig. 5. 
This code tree is obtained by writing along each branch the two digits of ^ 
representing the encoder output corresponding to the encoder input sequence 
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The non- systematic structure of the code allows for a 'quick- look- in' 
feature to be exploited (cf [1]). Quick-look decoding peimits recovery of 
the original information sequence from the hard-decisioned received data 
without resorting to sequential decoding methods simply by modulo-two 
addition of the received sequences. The quick-look decoding technique is 
illustrated in Fig. 6. Inverting the Q channel symbol at the encoder prior 
to transmission requires that the hard-decision made on the second symbol of 
every received channel pair be complemented (inverted) before the quick- 
look decoding can be performed. Since the taps for the two output symbols 

P and Q are identical, except for the second stage, then the original 
information sequence, D, c.Sin be reconstructed by modulo-two addition of 
the two symbol streams except for a one-bit delay. That is, the operation' 
has the form: 

Dj = Pj^3_ ® J = 0,1,2,... ,511 (3.1) 

This process will recover the complete information sequence except for the 
last bit received just before resetting the encoder. The s~ate of this 
bit cannot be predicted. However, the last bit of the original information 
sequence is the last bit of the frame synchronization word. When operating 
in the quick-look mode, it requires only that the last bit of the frame sync 
be ignored. The Reverse Playback Telemetry Program possesses the option 
of decoding a frame using the quick-look technique. Hard-decisions on the 
received sequence are obtained from the three-bit soft-decisions provided 
on the input digital tape. 

Here Q,^-, represents the logical compliment of 3,^.-, • Ideally this should 
equal in the absence of channel noise_^ We prefer the notation 

to emphasize that the quantities and in 2q- (3.1) are channel* 

outputs . 
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4,0 Frame Synchrotiizatlon Algtsrithm 

Due to the nature of sequential decoding^ it is necessary to synchronize 
the decoder with the input data before the decoding process can proceed. 
Synchronization is accomplished by locating the position of the sync word 
in each frame. 

Consider the problem of locating the sync vrord in n frames of data. 

It is highly probablff that the actual start bits of the data frames do 
not coincide exactly with the start bits of the physical frame records as 
they are stored on the digital data tape. This is demonstrated in Pig. 7 a. 

It is obvio.us that, in general, the n data frames will span n+1 tape frames. 

Let S. j denote the i'th symbol of the j'th tape frame (i = 0,1, ..*,1024; 
d 

j a 1,..., n+l). As noted previously, the coded symbol stream may be con- 
verted to the ox’iginal data stream by first inverting the hard-decision made 
on the second symbol of every channel pair and then perfonaing an exclxisive- 
or operation on the channel symbols. However, the correct pairing of the 
symbols for performing the ouldE-look decoding operation is not evident in 
the received data stream. This means that the quick-look decoding operation 
must be performed on the received data for the two possible combinations of 
channel symbol pairs, that is 


1. ) Sq^ J ® J ; Sg^ J © J ; . . . ; s^^gg^ ^ ® s^Qg^^ j 

2. ) ® Sg^j ; ® ® 


(4.1) 

(4.2) 


Now let Cj denote the i'th bit of the frame synchronization wrd (i =0,1,..., 
22). The recovered data streams, by So.ns. (4.1) and (4.2), are then each 
correlated with the sync word for each possible symc word location. The sync 
word location which gives the highest correlation is taken as the true sync word 
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location. This can be stated mathematically as follows: 

Define the correlation operator * according to the truth table of Pig. 
7b, This operation gives a value of -1 if the symbols to be correlated 
are different and +1 if the symbols ore the same. The operation thus per- 
formed is: 


max max ^ I * 

m«0,l; • • • > 511 k*0,l J«1 i *0 


\ j ^+k+2 I +1}, j 


(l*.3) 


Where, k represents the two possible channel symbol pairings and m represents 
the 512 possible locations of the start of the sync word. Note that the first 
subscript of each channel symbol term is enclosed in braces. This indicates 
the term is to be evaluated modulo 102b for the case k=0 and modulo 1025 


for the case k=l, i.e., 


0 ^n< 1023 
n > 102b 




1025 n-1025 


0 ^n< 102b 
n > 1025 


This takes into accoiint the possibility that the sync word may cross over 
a tape frame boundary as illustrated in Fig. 7c. It is easily see that if 
L bits of the sync word occur at the end of the tape record, then the re- 
maining 2b-L sync bits from the previous data frame will occur at the start 
of the tape record. The modulo operation allows the correlation to be per- 
formed on the L bits at the end of the record and the 2b-L bits at the be- 


ginning of the record. 
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After the maximum has been found, the values of h and m give the proper 

channel symbol pairing and location of the sync word starting point. The 

location of the start of a frame is easily deduced. 

4.1 Frame Synchronization Subroutine Description 

The subroutine STCNC performs the synchronization of the data for the 
main program. The call to the routine takes the form: 

CALL SYNC ( J, NBEC, IDTSK, IPBB, IWOHK, IHETA, iqU, NPRMS , ESNO, IFFR3; ISYCVr ) 
where the parameters are; 

J - is the first record to be used in locating the sync word. 

NEEC - the total number of records of input data. 

IDISK - is the logical unit number for the random access device (^6). 
IPRB - is a ( 512 , 9 ) array used for inputting the data from which 
the sync word is to be located. 

IWORK - is a (1025,9) array used as a work array. 

IBETA - is a (1024,8) array in which che unpacked sof^ -decisions 
are returned in proper synchronization. 

IQLA - is a (512,9) array used as a work array. - Upon subroutine 
completion it contains the quick-look outputs in proper 
frame sync . 

NFEMS - is the number of frames which are to be used in locating 
the sync word. (1 ± NFRMS £ 8) 

ESNO - is an eight word array in which the signal to noise ratio 
for each synchronized frame is returned. 
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IFFRT - is the number of frames synchronized and returned in 
IQLA and IRETA. (O <. IFFRT 1 8) 

ISYCW - is location of the frame sync word returned by 
the subroutine. 

The subroutine determines the location of the frame sync word in the 
received data sequence and then re-orders the sequence to obtain proper 
frame synchronization. Having acquired frame sync the subroutine computes 
the quich-look outputs. Purtheimore, due to the inverter on the Q channel 
symbol output of the encoder^ the complement of the Q received soft= 
decisions are computed by the SYNC subroutine providing the encoder outputs 
?, Q to the Fano subroutine in proper frame sync. 
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5.0 Fano Seouential Decoding Algorithm 

The sequential decoder accepts a sequence of demodulator output symbols, 

denoted by r * {r.}, and attempts to find a path through the code tree describ- 

<3 

ed in Fig. 5 , which has a high likelihood of having produced the sequence r. 

It does so by selecting a tentative path into the tree starting at the origin 
and following at each node the branch that best matches the appropriate segment 
of r. Whenever the path that the decoder is cturently following becomes too un- 
likely, a search is initiated for a better path. The attractive feature of a 
sequential decoder is its ability to determine the correct (transmitted) L-branch 
path while examining only a small fraction of the 2 (2^-1) branches that make up 
the votal set of 2^ possible paths. The actual number of branches ex ami ned in 
decoding a sequence is a random vairiable. If the channel noise is quite severe, 
the number of branches that must be examined (some more than once) by the decoder 
before accepting an L-bit path becomes very large. The examination of a branch, 
in either the forward or reverse direction, will be referred to as a computation. 

The data have been organized into frames of li88 information bits long. The 
32-stage shift register in the convolutional encoder is set to a prechosen value 
of all O's before transmitting each frame. After each frame is transmitted a pre- 
chosen sequence of 2^* bits is fed into the encoder. This, known sequence, called 
the tail sequence, serves to terminate the code tree kS channel symbols sifter the 
final node, thereby providing sufficient energy to reliably decode the last inform- 
ation bit. Furthermore, the known tail sequence provides a repeating pattern of 
channel symbols (independent of data) which is used for obtaining frame synchroni- 
zation on the received bit stream. Each telemetry frame is composed of 1024 channel 
symbols corresponding to 512 data bits [ 2 ] . 

The convolutional encoder output symbols are transmitted to the receiving 
station and are demodulated at the receiver output. Each output symbol represents 
the sum of the transmitted symbol and the system noise. Each demodulator output 
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is quantized to one of eight amplitude levels. These soft-decisions on the 
received sequence are made use of by the sequential decoder. The decoder 
contains a replica of the encoder in which it can try various possible trans- 
mitted sequences and compare its output with the channel output. The decoder 
attempts to find the correct transmitted path thro\:gh the code tree by making 
use of reliability information provided by the soft-decisions. In order to 
minimize the probability of decoding on the wrong path, the decoder must, in 
the decoding process, attempt to find the most probable path through the tree. 


By Bayes rule, the a posteriori probability of any path K nodes long is 

p(r|x^J^)p(x^J^) 


p(x^-^^lr) = 


p(r) 


( 5.1) 


.(J) 


where the vector r is the sequence of channel outputs, and x is the sequence 
of channel inputs corresponding to the jth path. Any monotonically increasing 
function of the a posteriori probability will suffice in characterizing the 
performance of the decoding process. Since the system noise is additive white 
Gaussian and the channel is memoryless ( 1 ) reduces to a product of probabilities, 
An equivalent measure of path likelihood is given by the cumulative path metric 
(up to level K). 

K I n r 1 j 

( 5.2) 


I 

i=l ‘ 


n 

I 

mal 




Here, r^^^; m=l,2,...,n is the m"**^ receiver output along the branch and xj^^ 
is the associated channel input sequence. The constant 3 is chosen to have a 
value such that, for the correct path, L, will increase on the average and, for 

Jo. 

incorrect paths, L, will decrease on the average. The Fano algorithm detects 
incorrect paths by comparing the cumulative path metric with a running thres- 
hold value. 
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Fig. 3 shows the relationship between L and the sequence of thresholds t 
separated by T^, as the decoding process evolves. As the decoder proceeds 
forward along a path, It chooses the most likely branch at each node and then 
tlght«.= the nmhlhS thre.hold eo thM It 1. never mre then heneath for 
the Gxirrent path. The flow diagram of Fig. 9 Illustrates the action of the 
Fano algorithm. If the path should fall below the threshold, the decoder goes 
in a backward search in which it searches all other paths in order of decreas- 
ing branch probability looking for one which does not fall below the present 
threshold T. If all those paths which originate above the threshold are ex- 
hausted, the threshold is reduced by Tq and the original path is retraced. If 
this path doss not remain above the sew threshold, the decoder Tearches back 
again looking for a path which does. This process is repeated, reducing the 
threshold by steps, until success is realized, at which point the decoder con- 
tinues forward into the decoding process. To eliminate the possibility of the 
decoder entering into a loop in the searchback mode the threshold la not allow- 
ed to increase when the decoder reaches previously explored branches. This 
action is insured by making use of a binary /ariable 6 which keeps track of when 
the running threshold may be tightened. The variable 9, initially 0, is set 
equal to 1 immediately following observation of a i”jnni'ng threshold violation 
on a forward look. As long as 9 remains equal to 1 the algorithm prevents 
tightening of T. As soon as the new node is examined (one never reached before) 
9 is reset to 0 and tightening ,ts again permitted [; 3], 

5.1 Branch Metric Computation 

For the Fano sequential decoder, the cumiilative path metric up to level K 

f \ ) 

in the code tree associated with the path x"' on the basis of the receiver out- 
put sequence r = (r^jrg,...) is given by (5 .2). The Branch Metric associated 
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with the branch along path is given by the inner summation of (l) 


-th 




m*l 


log« 




- B 


( 5 . 3 ) 


Where, 3 is a bias term chosen to insure that (l) increases in value when the 
decoder is following the correct path. The bias term will be taken as equal to 
the code rate R» 1/2. 

The Pono decoding algorithm makes use of soft-decisions on the channel 
outputs ( demodulator outputs ) . Each channel output is quantized to one of eight 
levels. To a first order approximation, the modulator-channel-demodulator 
cascade can be replaced by a constant, discrete, memoryless channel with two 
input letters, eight output letters, and transition probabilities{p(r^^^ 

Under the assumption of coherent BPSK modulation with matched filter reception 
appropriately normalized to the noise power the complete received signeJ. 


r(t)=s(t)+n(t) can be characterized os 


T SS X 

Im Im 



+ a 


im 


( 5 . 1 +) 


where is the transmitted binary symbol (.±1) and is a zero-mean, unit 


varian.ce Gaussiain random variable. It follows that 


while 




The channel transition urobabilities can now be calculated as follows: 




Im 


( 5 . 5 ) 

( 5 . 6 ) 


( 5 . 7 ) 


Hera , 




k=0,l,..., 8 are the input bin boundaries of a uniform, symmetric 


quantizer (Fig. 10) and r 


1 (k) 


2-m 


k=l,2,..., 8 are the output values of the 
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’ ■ ”t Po'V’P'^ta l = ta> 8 (5.8) 


quantizer. (Note: 2 q»-“, Eg»+<"). The quantity is given by (5.6) 

Hence, for quantized outputs (5. 6) becomes 

n-1 

Jo 

In probability «.* transmitting symbol x^. It is 

assumed that both binary symbols are equally likely to be channel inputs. 

Noting that there are two channel symbols per branch (n“2) explicit evalu- 
ation of the branch metric can now be performed, 

[ih ,1 

(5.9) 


4'"> 


2 

I 

m*l 


^(£^.4^’)- 1082 


log 






+ lOgg 


'i2 ^£2 
Pn(i‘»o) 


P(^ool^oo ) 


- 1 


(5.10) 


Evaluating each term of (5.10) 

JJ) 


log^ 






m=0 




log. 






;(k) 


for i=l,2, 

(5.11) 

,.■=11 


Since, can be one of eight possible values, denote to be the k quant- 

izer o\,tput value. Then, rewriting the right-hand side of (5.II) yields, 

8P(?i?’l4i’) 


log. 


'ai '^ii 


where represents tue complement of The branch metric expression now 


becomes 


= i°s. 


2p(f'jl>U<J>) 

■^•lOgg 

2p(4^’ 

^2^) 1 




(5.12) 
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Eo.n. (5.12)can be written more simply due to symmetry properties of the quantiser 
and resulting transition probabilities. 


Define: 


Branch Metric: 


TP(K) ^ 

TP(9-K) ^ 

TP(J) ^ p(r[^2)|jc(j)) 
TP(9-J) ^ 


‘12 


■ logg 


’ 2TP(K) 

tp(k)+tpT^kF 


+log. 


2 TP(J) 


TP(J)+TP(9-J) 


(5«13a) 

(5.13b) 

(5.13c) 

(5.13d) 

-1 

(5.1H) 


Eqn. (5.l4)holds for K»l,2,...,8; Js»l,2, . . . ,8. There are 6U possibilities arising 

from noting that individual branches are represented by a pair gf channel symbols 


each of which is quantised to one of eight possible values. Each value determined 
by(5.lU)is being pre-computed and entered into a table. The decoder mahea use 
of the table to look-up the branch metric value associated with a particular branch 
at the current node. The look-up table is structured as shown in Fig. 11 . 


5.2 Parameters Effecting Seouential Decoder Performance 

An important and useful parameter in characterising the performance of the 

Fano sequential decoding algorithm is the critical rate R- (often referred to as the 

u 

computational cutoff rate R _ in this context) which for coherent 3PSK modulation 

comp 

and equally probable signaling is given by (cf. [3 ], Chap. 6) 


RQ»l-logg[l+ 


/p{r|x*l}p{r|x=-l} dr] 


( 5.15) 


where p{r|x=i} is the aposteriori probability of having observed a receiver output 
value r during any one channel signaling interval given that the corresponding 
transmitted symbol was x=»i for i= +_ 1. The decoder makes use of soft decisions 


2k 


0 




on thfe receiver outputs, hence for quantized receiver outputs, the value of 
Rq must be computed according to 

g 

Rg * 1 - log 2 [l+ I /p{i^qj lx*l}p{f»qj lx»-l} ] (5.l6) 

j*l 

where now p{r»qj|x»i} is the a posteriori probability that the qxoantized 
receiver output assumes the J quantization level q^ ; J=l,2,...,8 during any 
signaling interval given that the corresponding transmitted symbol was x*i 
with i= + 1. and is given by (5-7) 

The convolutional code rate R, in terms of information bits per chaimel use, 
normalized to the critical rate Rg provides an indication of the behavior of the 
average nmber of computations required per decoded information digit. The para- 
meter H/Rg characterizes sequential decoder operation in terms of the channel 

properties and the quantizer characteristics. A nlot of R/R„ versus E /N- in dB 

0 so 

is shown in Fig. 12. For the system of interest R=l/2, and the quantizer is both 
symmetric and uniform with input bin boundary spacing A=0.5 (normalized in stand- 
ard deviations of the system noise). It is noted that for values of E^/Ng such 
that R/Rg<l the Fano algorithm will with high probability, decode a frame success- 
fully since the average number of computations required per decoded information 
bit is finite. Further, when E^/Ng is such that R/Rg^l the average number of 
compiitations required to decode a frame becomes unbounded. Table II contains 
selected values of E^/Ng in dB and the corresponding values of R/Rg are provided 
along with a tabulation of the channel transition probabilities for each case. 

A series of simulations were performed for the Fano decoder used in this 
program [ 9 ] • 
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5.3 Fane Algorithm Subroutine DeseriPtlon 

The Fano sequential decoding algorithm has been designed as a sub- 
routine which is cailed by the main driver program. The Fortran code 
appears in Appendix B. The program sequentially decodes frames of 102U 
channel symbols using soft-decisions on the demodulated channel output. 
Each frame corresponds to 512 information bits of which the last 2h bits 
is the known frame sync word. This is used by the decoder as a known 
tail.sequence . Decoding into the tail allows for reliable decoding of 
the last data information bit. The call to the routine takes the form- 
CALL FANO (ESNODB, IN, lOUT, ITCT, ACB, NSC, IDm,G, ITPPLG) 
Where, ESNODB « signal energy to noise density ratio in dB. 

IN ■ input array of 1024 soft-decisions that are input 
to the Fano decoder. 

lOUT = Output array of 512 decoded information bits, 

ITCT » on entrance to routine: maximum number of com- 
putations allowed per frame desired by tsser (defaiilt 
is 100,000) on exit from routine: number of for- 
ward and reverse moves made by the decoder in pro- 
cessing a frame. 

ACB = average number of computations per decoded bit . 

NSC = number of corrected channel symbols (ec.uals zero 
when frame is deleted) . 

IDFFLG = deleted frame flag: 

0 ■ frame was successfully decoded 

1 a frame deleted because the maximum number of computa- 

tions allowed per frame has been exceeded. 

2 “ invalid decoder parameter has been encountered, 

ITPFLG = flag indicating whether output is written to 
magnetic tape. 

The decoder calculates the channel transition probabilities from the 
signal-to-noise ratio read from the start of each telemetry frame. 
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The iinpliclt assumption Is that the received sequence has been quantized 
using a uniform quantizer, vherein the input bln boundaries have been 
normalized to the standard deviation of the noise power. The value of 
d*0,5 has been fo^^nd to be near optimum Having the transition pro- 

babilities the decoder proceeds to calculate the branch metrics according 
to f5.ll»). The branch metrics are scaled to allow integer computations 
in the decoding process. The branch metric values are stored in two arrays 
as shown in Pig. 11. These arrays are accessed by an index which is dete- 
mined from the pair of received channel symbols representing a particu- 
lar branch in the code tre-s. The IN array contains 102U values in the 
range 1 to 8, representing the quantizer levels, the received channel 
symbols have been assigned. 

At the start of each frame the decoder is initially set to zero. 

The decoding process continues until either the frame is decoded success- 
fully or the frame is deleted. Upon successfully decoding a frame the 
subroutine returns the average number of computations per bit. This is 
calculated by dividing the number of forward and reverse moves by the 
number of information bits in a frame. Also, the nusber of corrected 
cbann,el symbols is returned by the subroutine . This is determined by 
comparing the branch labels determined by the replica of the encoder, 
which is embedded in the Fano algorithm, with hard-decisions made on 
the channel output symbols contained in the IN array. A frame is deleted 
if the maximtaa number of computations is exceeded in attenroting to de- 
code a particular frame. If the Printout option is specified, the bits 
declared to be decoded will be printed out . Undecoded bits will be 
indicated by a "9" acting as a place holder in the data bit sequence. 

The frame sync word (bits U89-512) will be printed out for reference. 

For example, when a frame deletion is encountered the partially decoded 
frame will appear on the print-out as follows - 


5 


FRAHE WAS DECODED USING FANO ALGORITHM* 

FRAME WAS DELETED. 

NUMBER OF SYMBOL ERRORS CORRECTED® 0 

AVERAGE NUMBER OF COMPUTATIONS PER BIT® 

<* BITS DECODED/MAX. ♦ COMP. PER FRAME) 

DECODED DATA- 

1001100100110111000000000000000000000000000000000000000000000000 

0000000000000000000000000000000000000000000000000000000000000000 

0000000000000000000000000000000000000000000000000000000000000000 

0000000000000000000000000000000000000000000000001000000010000001 

0101010000000000000000000000000000000000000000000000000000000000 

0000000000000000000000000000000000000000000000000000000000000000 

0000000000000000000000000000000000000000000000000000000000000101 

0011001001010001010100001199999999999999999999999999999999999999 

TAIL SEQUENCE SHOULD BE- 111110001100010101001001 

The purpose of printing out deleted frames is to allow the user 
to re-atteapt decoding the particular frame after increasing the 
maximum number of computations allowed per frame. 3y making several 
decoding trials on a difficult frame the user may be able to deter- 
mine which bits in the partially decoded frame may possess some degree 
of reliability as the decoder progresses into the known tail sequence. 

It should be noted, of course, that the accuracy of any and all bits 
in a partially decoded frame is highly questionable. Furthermore, 
computer CPU processing time will increase considerably on an IBM 
360/67 as the number of computations parameter is increased for frames 
containing a significant number of channel symbol errors. If the Tape 
Ctitput option is specified^ a deleted frame will be recorded on tape 
such that bits 193-680 are all zeros and bits 68I-70U of a Data Block 
contain the frame sync word. The Data Block is a data record contained 
in the Intermediate Data Record (IDB) file. Frame deletion will also 
occur if the decoder encounters an invalid parameter during execution, 
and in the print-out made an error message to that effect is printed. 
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LOGICAL UNIT NUMBER 

INPUT-OUTPUT DEVICE 

1 

user terminal input 

2 

user terminal output 

3 

print output 

^ n 

tape input 

5 

tape output 

6 

random access input and 
output 


Table 

Logical I/O Device 


1 

Assignments 
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TRANSITION PROBABILITIES = 0.53457 0.18675 0.14011 0.08229' 0.03783 0.01361 0.00383 0.00101 




PIONEER VENUS MULTIPROBE MISSION 

TELEMETRY DATA RECOVERY AND INTERFEROMETRY EXPERIMENT 
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Flowchart for Main Prograr: 
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RUN -LOAD+ANOEJFILEl l=*SOURCE* 2«>|CSINK5|: 3«^cPRINTii< 5»*T0:« 6»-D 

^EXECUTION BEGINS 

IS DATA RECORDED ON TAPE IN FORWARD CW? REVERSE DIRECTION (FWD OR REV)? 
REV 

DO YOU WANT TAPE OUTPUT OR PRINT OUTPUT (TAPE OR PRNT)? 

PRNT 

DO KOU WANT PRINTOUT TO THE TERMINAL OR TO THE LINE PRINTER 
(TERM OR LPRNT)? 

TERM 

WHICH PHASE REFERENCE DO YOU WANT- DATA OR DATA-BAR (D OR DB)? 

DB 

IS THE DATA TO BE DECODED OR IS RAW SYMBOL OUTPUT DESIRED (DEC OR RAW/'i 
RAW 

DO YOU WANT RAW SYMBOLS TO BE SYNCHRONIZED (1024) IN A DATA FRAME OR 

THOSE (102S) CORRESPONDING TO A TAPE RECORD (SYNC OR TREC)7 

TREC 

DO YOU WANT ALL THE FRAMES PRINTED OUT (Y OR N)? 

N 

DO YOU WANT TO SPECIFY THE FRAMES TO BE PROCESSED ACCORDING TO 
TIME TAG OR FRAME NUMBER (TIME OR NUMB)? 

NUMB 


SPECIFY BY NUMBER THE FRAMES WHICH ARE TO BE PROCESSED. 
THE LAST DIGIT OF EACH NUMBER MUST END UP IN COLUMN 5. 
INDICATE YOU ARE DONE BY HITTING A RETURN. 

20 

21 


ENTER DSS NUMBER (14 OR 43)- 
14 


ENTER PROBE ID (SPlfSP2fSP3 OR LP)- 
SPl 
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278 RECORDS HAVE BEEN READ FROM INPUT TAPE 


PIONEER VENUS REVERSE PLAYBACK TELEMETRY PROGRAM 
PROBE IDENTIFICATION- SPl 
DSS NUMBER- DSS 14 
PLAYBACK DAY OF YEAR- 89 


3U 


Figure 3 

Sanple Outpu- of Program 




DATA START TIME- 
HOURS- 17 
MINUTES- 53 
SECONDS- 39,680 


1 

ov 


^ I k 
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PRINTOUT OF RAW CHANNEL SYMBOLS FOR TAPE FRAME (RECORD) 20 

DATA START TIME- 
HOURS- 17 
MINUTES- 53 
SECONDS- 1,680 


DATA STOP TIME- 
HOURS- 17 
MINUTES- 52 
SECONDS- 59,680 

SIGNAL-TO-NOISE RATIO <DB)= 6,703125000 


RAW CHANNEL SYMBOLS (SOFT-DECISIONS IN OCTAL FORMAT)- 


3524352525352434243535240426251434242424341425363524261635353525 

1626351526150534061517263425252504001425142624251514252634204235 

4541264463323432663242462142524415341435315240542712434434204434 

3255016245243516141517352415241534262435352405261615342534342525 

3505172525251514262425141535261415352414353425253435241000062405 

1425242425350624362574063525263624261634353616142626262535163414 

3434162524151424243426152434242616361526242426251525350416243526 

1514252414342626342426743414142514162424242424350425242434253625 

2435251434371514253524343434152415353414242634343535351625143535 

2534753514153635143534253532511455633322434253344264353444211344 

4344122433515463424466434343006003312313524325325616263145352434 


3346324420552254326423552731553526235622551434343524353525253635 

1626252634351535242424343424251424353625342435252515352525353424 

2534253724151425160434261524343626143425351534242414242426270626 

1624242634352336241534341426242636042424242526142424142425251524 

15352605051615261424241424242536041524352624163534342525350505141 


r:' 


RAW CHANNEL SYMBOLS (HARD-DECISIONS IN BINARY FORMAT)- 

1010101010101010101010101010101010101010101010101010101010101010 
1010101010101010101010101010101010111010101010101010101010110110 
0001100001111011001101001101010010101010110101001011010010110010 
1100110100101010101010101010101010101010101010101010101010101010 
1010101010101010101010101010101010101010101010101010101111101010 
1010101010101010101000101010101010101010101010101010101010101010 
1010101010101010101010101010101010101010101010101010101010101010 
1010101010101010101010001010101010101010101010101010101010101010 
1010101010101010101010101010101010101010101010101010101010101010 
1010001010101010101010101011011000011111010101100100101.000111100 
0100111011010001010000010101110111111111010110110010101100101010 
1100110011001100110011001011001010110011001010101010101010101010 
1010101010101010101010101010101010101010101010101010101010101010 
1010101010101010101010101010101010101010101010101010101010101010 
1010101010101010101010101010101010101010101010101010101010101010 
loioioioioioioici 01 0101010101010101010101 0101 0101010101 01 01010101 
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Figure 3 - cont 
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PRINTOUT OF RAW CHANNEL SYMBOLS FOR TA^'E FRAME (RECORD) 

DATA START TIME- 
HOURS- 17 
MINUTES- 52 
SECONDS- 59.680 


DATA STOP TIME- 
HOURS- 17 
MINUTES- 52 
SECONDS- 57.680 

SIGNAL-TO-NOISE RATIO (DB)=» 6.703125000 


RAW CHANNEL SYMBOLS (SOFT-DECISIONS IN OCTAL »=’ORMAT)'- 

1634252425241535340536141424152436252525353434243534372526243634 

2525352534142505350634343425172435241535143435253525251504334657 

36436335216542466223541435263525134033S356262351426252424U352713 

4244052243425354342415161315241525241435353536363525143525343635 

3436262435351405242504142435363416362526242425243614361424252534 

1535353624353634242534263435253435153434353534372524352714241435 

1424262526243535243526242424343537143434253415352434353506253524 

1516263414153535361670342525252534342424353425241524262535262434 

2434241514362535363634353634263425151434142515162524243414243525 

3535J454 152363243427362424224015476 1232363524134534624254533135^ 

42J51224236I24442525545524353235512333206165622462612424354225343 

54242521636^2242653433660435235425251365150534143434063735353625 

2535353625351405272535101425243424342535052626341436242504353415 

2524253625243634261515242514342425253434353634353515253425351527 

2424343515240534362636253525152725261535343424242426241536341735 

24152424243615250524252536151614152434362515163434361414162636252 


RAW CHANNEL SYMBOLS (HARD-DECISIONS IN BINARY FORMAT)- 

10101010 10 10101 OlOlOlOlOlOlOlOlOlOlOlOlOlOiOlOlOlOlOlOlOlOlOlOlO 
1010101010101010101010101010101010101010101010101010101010110000 
1001011011000100011100101010101011011101001011010101010100101011 
0100101101010100101010101010101010101010101010101010101010101010 
1010101010101010101010101010101010101010101010101010101010101010 
10 10101010101010101010101010101010101010101010101010101010101010 
1010101010101010101010101010101010101010101010101010101010101010 
1010101010101010101001101010101010101010101010101010101010101010 
1010101010101010101010101010101010101010101010101010101010101010 
1010101010101010101010101011011000011111010101100100101000111100 
0100111011010001010000010101110011111110100011001011010100110101 
0010101 101001 10100101 10010101 iOOlOlO 11 00101010101010101010101010 
10101010101010101010101 110101010101 0101 0101010101010101010101010 
1010101010101010101010101010101010101010101010101010101010101010 
1010101010101010101010101010101010101010101010101010101010101010 
10101010101010101010101010101010101010101010101010101010101010101 


PROGRAM SUMMARY 

TOTAL NUMBER OF FRAMES PROCESSED= 
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Figurs 3 - cont. 


c 


5 

: . ■) 

. . Olr . ... 

ARE YOU DONE (Y OR N)T 
M 

DO YOU WANT TAPE OUTPUT OR PRINT OUTPUT (TAPE OR PRNT)7 
PRNT 

DO YOU WANT PRINTOUT TO THE TERMINAL OR TO THE LINE PRINTER 
(TERM OR LPRNT)? 

TERM 


WHICH PHASE REFERENCE DO YOU WANT- DATA OR DATA-SAR (D OR DB)7 
DB 

IS THE DATA TO BE DECODED OR IS RAW SYMBOL OUTPUT DESIRED (DEC OR RAW)? 
RAW 

DO YOU WANT RAW SYMBOLS TO BE SYNCHRONIZED (1024) IN A DATA FRAME OR 
THOSE (1025) CORRESPONDING TO A TAPE RECORD (SYNC OR TREO? 

SYNC 

DO YOU WANT ALL THE FRAMES PRINTED OUT (Y OR N)? 

N 

DO YOU WANT TO SPECIFY THE FRAMES TO BE PROCESSED ACCORDING TO 
TIME TAG OR FRAME NUMBER (TIME OR NUMB)? 

NUMB 

SPECIFY BY NUMBER THE FRAMES WHICH ARE TO BE PROCESSED. 

THE LAST DIGIT OF EACH NUMBER MUST END UP IN COLUMN 5. 

INDICATE YOU ARE DONE BY HITTING A RETURN. 

20 


PIONEER UENUS REVERSE PLAYBACK TELEMETRY PROGRAM 


PROBE IDENTIFICATION- SPl 


DSS NUMBER- DSS 14 


PLAYBACK DAY OF YEAR- 89 


DATA START TIME- 
HOURS- 17 
MINUTES- 53 
SECONDS- 39.680 


PRINTOUT OF RAW CHANNEL SYMBOLS FOR SYNCHRONIZED DATA FRAME 20 

DATA START TIME- 
HOURS- 17 
MINUTES- 53 
SECONDS- 1,680 


DATA STOP TIME- 
HOURS- 17 
MINUTES- 52 
SECONDS- 59.680 
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Figure 3 - cont. 


SYNC WORD LOCATION- 326 
SIGNAL-TO-NOISE RATIO <DB)» 6.703125000 






SYNCHRONIZED RAW CHANNEL SYMBOLS (SOFT-DECISIONS IN OCTAL FORMAT)- 

1546342446643434300600331231352432532561626314535243433463244205 

5225432642355273155352623562255143434352435352525363516262526343 

5153524242434342425142435362534243525251535252535342425342537241 

5142516043426152434362614342535153424241424242627062616242426343 

5253624153434142624263604242424252614242414242525152415352605051 

6152614242414242425360415243526241635343425253505051416342524252 

4153534053614142415243625252535343424353437252624363425253525341 

4250535063434342517243524153514343525352525150433465736436335216 

5424662235414352635251340335356262351426252424535271342440522434 

2535434241516151524152524143535353636352514352534363534362624353 

5140524250414243536341636252624242524361436142425253415353536243 

5363424253426343525343515343435353437252435271424143514242625262 

4353524352624242434353714343425341535243435350625352415162634141 

5353536167034252525253434242435342524152426253526243424342415143 

6253536363435363426342515143414251516252424341424352535353434152 

5352434273624242240154761232363524134534624254533135542451224236 


SYNCHRONIZED RAW CHANNEL SYMBOLS (HARD-DECISIONS IN BINARY FORMAT)- 

1000101000001010111011111111101011011001010110010101011001100110 
011001100110010110010101100110010.1010101010101010101010101010101 
0101010101010101010101010101010101010101010101010101010101010101 
0101010101010101010101010101010101010101010101010101010101010101 
0101010101010101010101010101010101010101010101010101010101010101 
0101010101010101010101010101010101010101010101010101010101010101 
0101010101010101010101010101010101010101010101010101010101010101 
OlOlOlOlOlOlOlOlOlOlOlOlOlOlOiOlOlOlOlOlOlOlOlOl 10000100101 101 10 
0010001110010101010101101110100101101010101010010101101001011010 
1010010101010101010101010101010101010101010101010101010101010101 
OlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOiOlOlOl 0101 0101010101010101010101 
0101010101010101010101010101010101010101010101010101010101010101 
OlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOlOl 
0101010100110101010101010101010101010101010101010101010101010101 
0101010101010101010101010101010101010101010101010101010101010101 
0101010101010101101100001111101010110010010100011110001001110110 


PROGRAM SUMMARY 

total number of frames processed* 1 

ARE YOU done (Y OR N>? 

N 

DO YOU WANT TAPE OUTPUT OR PRINT OUTPUT (TAPE OR PRNT)? 

PRNT 

DO YOU WANT PRINTOUT TO THE TERMINAL OR TO THE LINE PRINTER 
(TERM OR LPRNT)? 

TERM 

WHICH PHASE reference DO YOU WANT- DATA OR DATA-BAR (D OR DB)7 
DB 
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Figure 3 -cont 


OF 


!. rr'r^ b 

hCO i Q<J^'‘-5rY 


IS THE DATA TO BE DECODED OR IS RAW SYMBOL OUTPUT DESIRED (DEC OR RAW)? 
DEC 

DO YOU WANT THE FAMO DECODER USED OR THE OUICK-LOOK OUTPUTS 
(FANO OR QD? 

FANO 

INPUT MAXIMUM NUMBER OF COMPUTATIONS PER FRAME. 

THE LAST DIGIT OF THE NUMBER MUST END UP IN COLUMN 8. 

0 INDICATES DEFAULT»100f000 COMPUTATIONS. 

125000 

DO YOU WANT ALL THE FRAMES PRINTED OUT (Y OR N)? 

N 

DO YOU WANT TO SPECIFY THE FRAMES TO BE PROCESSED ACCORDING TO 
TIME TAG OR FRAME NUMBER (TIME OR NUMB)? 

NUV 

DO YOU WANT TO SPECIFY THE FRAMES TO BE PROCESSED ACCORDING TO 
TIME TAG OR FRAME NUMBER (TIME OR NUMB)? 

NUMB 

SPECIFY BY NUMBER THE FRAMES WHICH ARE TO BE PROCESSED. 

THE LAST DIGIT OF EACH NUMBER MUST END UP IN COLUMN 5. 

INDICATE YOU ARE DONE BY HITTING A RETURN. 
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PIONEER UENUS REVERSE PLAYBACK TELEMETRY PROGRAM 


PROBE IDENTIFICATION- SPl 


DSS NUMBER- DSS 14 


PLAYBACK DAY OF YEAR- 89 


DATA START TIME- 
HOURS- 17 
MINUTES- 53 
SECONDS- 39.680 


PRINTOUT OF DECODED DATA FOR FRAME 277 

DATA START TIME- 
HOURS- 17 
MINUTES- 44 
SECONDS- 27.677 
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Figure 3 - cont. 



r, . \ ^ 

yt , \ b. ^ If ’ ^ 

OF FOOK QUALSTV: 


DATA STOP TIME- 
HOURS- 17 
MINUTES- 44 
SECONDS- 25,677 

SYNC WORD LOCATION* 326 

SIGNAL-TO-NOISE RATIO <DB)« 0,5000000000 

FRAME WAS DECODED USING FANG ALGORITHM, 

FRAME WAS DECODED SUCCESSFULLY, 

NUMBER OF SYMBOL ERRORS CORRECTED* 1 

AUERAGE NUMBER OF COMPUTATIONS PER BIT* 1,0000000 

DECODED DATA- 

1001 1001001 1 1010000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000011111111 
0000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000111110001100010101001001 


PROGRAM SUMMARY 

TOTAL NUMBER OF FRAMES PROCESSED* 1 

TOTAL NUMBER OF FRAMES DELETED* 0 

DELETION RATE* 0,0 

SYMBOL ERROR RATE* 0,000977 

ARE YOU DONE (Y OR N)7 
N 

DO YOU WANT TAPE OUTPUT OR PRINT OUTPUT (TAPE OR PRNT)? 

TAPE 

DO YOU WANT PRINTOUT TO THE TERMINAL OR TO THE LINE PRINTER 
(TERM OR LPRNT)? 

TERM 

WHICH PHASE REFERENCE DO YOU WANT- DATA OR DATA-BAR (D OR DB)? 
DB 

DO YOU WANT THE FANO DECODER USED OR THE OUICK-LOOK OUTPUTS 
(FANO OR QD? 

FANO 

INPUT MAXIMUM NUMBER OF COMPUTATIONS PER FRAME. 

THE LAST DIGIT OF THE NUMBER MUST END UP IN COLUMN S. 

0 INDICATES DEFAULT=100f0O0 COMPUTATIONS. 


ENTER TAPE SEQUENCE NUMBER- 
12 


1.0 


rigTire 3 - oont. 


OF PO^-* 


ENTER YEAR DATA WAS RECORDED <E.Q. 78)- 
78 


PIONEER yENUS REVERSE PLAYBACK TELEMETRY PROGRAM 


PROBE IDENTIFICATION- SPl 


DSS NUMBER- DSS 14 


PLAYBACK DAY OF YEAR- 89 


DATA START TIME- 
HOURS- 17 
MINUTES- 47 
SECONDS- 7.678 


FRAMES DECODED USING FANO ALGORITHM. 

SYNC WORD LOCATION* 326 
SIGNAL-TO-NOISE RATIO <DB>= 6.968750000 

SIGNAL-TO-NOISE RATIO (DB>= 6.890625000 

SIGNAL-TO-NOISE RATIO <DB)= 6.890625000 

SIGNAL-TO-NOISE RATIO (DB>= 6.890625000 

SIGNAL-TO-NOISE RATIO (DB>= 6.890625000 

SIGNAL-TO-NOISE RATIO <DB)= 0.5000000000 
SIGNAL-TO-NOISE RATIO (DB)= 0.5000000000 
SIGNAL-TO-NOISE RATIO <DB)= 0.5000000000 
FRAMES DECODED USING FANO ALGORITHM. 

SYNC WORD LOCATION* 326 

SIGNAL-TO-NOISE RATIO <DB)= 0.5000000000 


PROGRAM SUMMARY 

TOTAL NUMBER OF FRAMES PROCESSED* 9 

TOTAL NUMBER OF FRAMES DELETED* 0 

DELETION RATE* 0.0 

SYMBOL ERROR RATE* 0.001411 


ARE YOU DONE (Y OR N)? 




Y 

^EXECUTION TERMINATED 


FigTire 3 - cont. 


EXCLUSIVE 


OF POOH 






































DATA niAKE n 


DATA FRAl^ n***! 


TAPE RECOBD n 


TAPE RECOBD n+1 


a. ) Overlap of data frames on tape records . 


X 

Y 

(X»Y) 

0 

0 

1 

0 

1 

-1 

■ 

0 

-1 

■ 

1 

1 


b.) Truth table for correlation operator. 


1023 


DATA FRAME 
n-1 


I SYRC WORD 
I FOR FRAME 

-L. ■ — nrJ. 


DATA FRAME 
n 


SYNC WORD 

FOR FRAME 
n 
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c.) Overlap of sync vord on tape record boundary. 


Figure 7 

Data Frame and Tape Record Relationship 










at Node I Threshold 


Action ( X indicates *'sct d m i 




0 0 

look at 1 

point to 1 


0 0 

look at 2 

point to 2 

set r* - 

-A 0 

look at 2 

point to 3 


-A 0 

look at 4 

X look at 2 

point to 2 

—A I 

look at 5 

X look at 1 

set r** 0 

0 1 

look at 3 

point to 3 


0 t 

look at 4 

point to 4 

set 0 M 0 

0 0 

look at 6 

X look at 3 

point to 3 

0 t 

look at 7 

point to 7 

set 9 M 0 

0 0 

look at 8 

X look at 3 

point to 3 

0 1 

look at 2 

point to 2 


0 1 

took at 5 

point to 5 

set 9 » 0 

0 0 

look at 9 

X look at 2 

point to 2 

0 1 

look at 1 

point to 1 


0 1 

look at 10 

point to 10 

set 9 M 0 

0 0 

took at 1 1 

point to 1 1 

set r- - 

-A 1 0 

look at 12 

point 10 12 



Figure 

a 


An exaiaple of "biie seGjrciiiag procedure 

used iDy 

the Fane 

Algorithm 






0 



Figure 9 

Flow diagram for the Fano Algori tzim 
employing a flag state for threshold tightening 














Qp pv'W** 





SCALED 

BRANCH 

METRIC 

VALUES MOO 

ai 



988 

949 

824 

409 

-589 

-2173 

-4046 

-7040 


.949 

911 

785 

370 

-628 

-2212 

-4085 

-7079 

0 

824 

785 

659 

244 

-754 

-2338 

-4211 

-7204 


409 

370 

244 

-170 

-1169 

-2753 

-4626 

-7619 

Symbol 1 

-589 


“-7S4 — 


-2168 

-S7S2 

-562s 

-B618 

-2173 

-2212 

-2338 

-2753 

-3752 

-5336 

-7209 

-10202 


-4046 

-4085 

-4211 

-4626 

-5625 

-7209 

-9082 

-12075 

X 

-7040 

-7079 

-7204 

-7619 

-8618 

-10202 - 

■12075 

-15069 



0 


Symbol 2 

1 







SCALED 

BRANCH METRIC 

VALUES MOl 


• 


-7040 

-4046 

-2173 

-589 

409 

824 

949 

988 


-7079 

-4085 

-2212 

-628 

370 

785 • 

911 

949 

0 

-7204 

-4211 

-2338 

-754 

244 

659 

785 

824 


. -7619 

-4626 

-2753 

-1169 

-170 

244 

370 

409 

Symbol 1 

=^618 


-3752 

-2168 

-1169 

^733 

“628 

— 

-10202 

“7209 

-5336 

-3752 

-2753 

-2338 

-2212 

-2173 


-12075 

-9082 

-7209 

-5625 

-4626 

-42,11 

-4085 

-4046 

1 

-15069 

-12075 

-10202 

-6618 

-7619 

-7204 

-7079 

-7040 



0 


Symbol 2 

1 




Figure n 

Branch Metric Array Table 


— 

--q.QO 


-E.oa a.QO 2.aa 

E^/Kq, in dB 


H. QQ 


6.00 


0.00 


SequentieJ Decoding Parameters for Octal Channel 
Figure 12 
50 
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Appendix A 

Main Driver Program Fortran Code 


OF POOR QUALSTY 


I 


1 C 


3 

A 


I 


I 

1 

I 

I 


6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


22 

23 

24 


05 

Aim 


26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 


C THIS PROGRAM READS A DIGITAL TAPE WHICH CONTAINS THE 
C SOFT-DECISIONS FROM THE DEEP SPACE NETWORK'S SYMBOL 
C SYNCHRONIZER FOR THE DATA WHICH WILL BE RECORDED 
C DURING THE OFF-LINE PROCESSING OF THE PIONEER- 
C VENUS PRE-DETECTION RECORDINGS. ALSO INCLUDED ON 
C THE TAPE IS TIMING AND IDENTIFICATION INFORMATION. THE 
C PROGRAM FINDS THE 24 BIT SYNC WORD OF EACH TELEMETRY 
C FRAME ON THE TAPE» THEN SEQUENTIALLY DECODES EACH FRAMED 
C AND WRITES A DIGITAL TAPE WHICH WILL CONTAIN THE DECODED 
C DATA ALONG WITH TIMING INFORMATION AND AN IDENTIFICATION 
C RECORD. 

C 

C SPECIFY THE VARIABLES. DOUBLE PRECISION VARIABLES 
C ARE USED FOR HOLDING CHARACTERS ONLY. 

C 

DOUBLE PRECISION DSS14r DSS43f DSSrOl >02»0UT0PT 
C 

C HALFWORD INTEGERS WILL BE USED FOR THE MOST PART 
C IN PROCESSING THE DATA. 

C 

COMMON IPHASE 

INTEQER*2 lASC t IND f ISPA f ILPA f ISEQ r IFR ? IFRl ? IFRB ? 
*IWORK r IRETA f IQLA r IFLBR » ILMA f I IDA » IF . IWD f lEV f I YEAR f 
*I YEARS f IHALF r IDATAf IDAYHr IDAYT r IDAYU r IPBID > 

*IFILL t M0R2 f MAND2 » MLSL2 r IHEXl 
C 

C DIMENSION ARRAYS USED BY PROGRAM. 

C 

DIMENSION IFRAMEdOO) fSFRAME(100f3) f IFR<521) f 
NcIFRKSZDf 

:<«IFRB<521>9) f IW0RK(1025f 9) rIRETA<1024f8) f IQLA<512f9) » 
*ESN0<8) »IFLBR<608) »IHALF(2) >NSC<8) fACB<8> . 

*IDFFLG<8) >ITCT<8) 

C 

C EQUIVALENCE VARIABLES. THIS IS DONE FOR THE MOST PART 


37 C TO FACILITATE LOGICAL OPERATIONS AND READING AND 

38 C WRITING FROM TAPE. 

39 C 

40 EQUIVALENCE < IVALr IFR <3 ) ) r < IFULL» IHALF < 1 ) ) f < IFLBR ( 7 ) » 


41 RIVALS) 

42 C 

43 C SET UP DATA CONSTANTS USED IN PROGRAM. FOR THE MOST PART 

44 C THESE ARE ALL ALPHANUMERIC CONSTANTS. 

45 DATA lTWDBfIREVB»ITP»IPR»IDECAf IRAWAfIFANAflQLAAf 


46 *IYES » INOf ISPA » ILPA f ILPIA » ISPIA » ISP2A r ISP3A f 

47 !«DSS14»DSS43»ILMAfIIDA.IWDf lEVflYEARf IFILLrITMA. 

48 *INMBAf IHEXlf I8Af I9Af IlOAf IllA/'FWD '.'REV '. 

49 *'TAPE' f 'PRNT' . 'DEC '.'RAW '.'FANO'.'QL '. 

50 *'Y '.'N ' . 'SP' . 'LP' . 'LP '.'SPl '.'SP2 '. 

51 :«'SP3 '.'DSS 14 '.'DSS 43 ' . Z4C4D . Z4944 . Z6744 . 

52 *24556. '78 ' . Z8888 . ' TIME ' . 'NUMB' . ZFFFF . Z00003800 . 

53 *200003900. Z00313000.Z00313100/ 

54 DATA IDDD.IDBAR.01.02/'D '.'DB '.'TERM '.'LPRNT '/ 

55 DATA IRTF.IRDF/'TREC' . 'SYNC'/ 

56 C 

57 C FORMAT STATEMENTS USED BY PROGRAM 

58 C 


OF poos QUAUfY 


59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 


C 

C 

C 

C 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


10 FORMAT <lHOf 'IS DATA RECORDED ON TAPE IN FORWARD S 
)*c"OR REVERSE DIRECTION (FWD OR REV)?') 

11 FORMAT <lHOf 'IS THE DATA TO BE DECODED OR IS RAW 'r 
*' SYMBOL OUTPUT DESIRED (DEC OR RAW)?') 

12 FORMAT <lHOr' DO YOU WANT THE FANO DECODER USED 0R'» 
t' THE QUICK-LOOK OUTPUTS'! 

*/lXr'(FANO OR QL)?') 

13 FORMAT (IHOf 'DO YOU WANT TAPE OUTPUT OR PRINT 'f 
*' OUTPUT (TAPE OR PRNT)?') 

14 FORMAT (IHOf 'SPECIFY BY NUMBER THE FRAMES WHICH ARE TO '» 
*'BE PROCESSED* '» 

*/1Xf'THE last digit of EACH NUMBER MUST END 'r 
:tc'UP IN COLUMN 5*'f 

*/1Xf' INDICATE YOU ARE DONE BY HITTING A RETURN.') 

15 FORMAT (IHOf' DO YOU WANT ALL THE FRAMES PRINTED 'r 
t'OUT (Y OR N)?') 

16 FORMATdHOf 'DO YOU WANT TO SPECIFY THE FRAMES TO 'r 
*'BE PROCESSED ACCORDING TO 't 

*/lXr'TIME TAG OR FRAME NUMBER (TIME OR NUMB)?') 

17 FORMAT (IHOf 'ENTER DSS NUMBER (14 OR 43)-') 

18 FORMAT( IHOf 'ENTER PROBE ID (SP1fSP2fSP3 OR LP)-') 

19 FORMAT( IHOf 'ENTER TAPE SEQUENCE NUMBER-') 

20 FORMAT (A4) 

21 FORMAT (A2) 

22 FORMAT (15) 

23 FORMAT( IHOf 'HOURS- ') 

24 FORMAT (IHOf 'MINUTES- ') 

25 FORMAT (IHOf 'SECONDS- ') 

26 F0RMAT(F12.5) 

27 FORMAT (A2f ID 

28 FORMAT (4(128A2)f9A2) 

29 FORMAT (1H1f5Xf 'PIONEER VENUS REVERSE PLAYBACK TELEMETRY 
4c' PROGRAM' f/) 


30 FORMAT (IHOfSXf 'PROBE IDENTIFICATION- 

31 FORMAT (1H0f5Xf 'DSS NUMBER- 'fA8f/) 


' fA4f/) 


32 FORMAT( IHOfSXf 'DATA START TIME- ' f /lOX f ' HOURS- 
*/10Xf 'MINUTES- ' f I2f/10Xf 'SECONDS- 'fF6.3//) 


I2i 




I 


I 


I 


1 


I 


I 


I 


I 

I 


] 


\ 

i 


t 


I 


119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 


Ui %■ L ,. .I'l k? 

Of' ■ ' 

33 FORMAT (IHOr 'ENTER YEAR DATA WAS RECORDED <E»Q. 78)-'; 

C 

' 34 FORMAT (lH0r/20Xf 'PROGRAM SUMMARY' r//6X» 

*' TOTAL NUMBER 'f 

*'0F FRAMES PROCESSED* ' r IS r/6Xr 'TOTAL 'f 
NUMBER OF FRAMES 'f 

♦ 'DELETED*- ' f I5f/6Xf 'DELETION RATE- ' fFIO. 6f 

♦/6Xf 'SYMBOL ERROR RATE-' fF10»6) 

C 

35 FORMAT (1H0f/22Xf' PROGRAM SUMMARY' f//6Xf 'TOTAL NUMBER 'f 
♦'OF FRAMES PROCESSED* 'fIS) 

C 

36 FORMAT (IHOf 'ARE YOU DONE <Y OR N)7') 

C 

37 F0RMAT<4<128A2)f96A2) 

C 

38 FORMAT< IHOf 'INPUT MAXIMUM NUMBER OF COMPUTATIONS' f 

PER FRAME* 

♦/IXf'THE LAST DIGIT OF THE NUMBER MUST'f 
♦' END UP IN COLUMN 8«'f 

♦/IXf'O indicates default- 1 00 f 000 COMPUTATIONS.') 

c 

39 FORMAT (18) 

C 

40 FORMAT (12) 

C 

41 FORMAT (IHOf 'ERROR OCCURRED WHEN TAPE WAS READ.'f 
♦/IXf'TAPE is incorrectly formatted or tape DATASET' f 

♦/IXf'WAS incorrectly specified.') 

c 

42 FORMAT (IHOf 'RANDOM ACCESS FILE COULD NOT BE OPENED. 'f 
♦/IXf'THE PROBLEM IS MOST LIKELY DUE TO INSUFFICIENT TEMPORARY' f 
♦/IXf 'STORAGE ALLOCATION OR INCORRECT DATASET DEFINITION.') 

C 

43 FORMAT( IHOf 'WHICH PHASE REFERENCE DO YOU WANT- 'f 
♦'DATA OR DATA-BAR (D OR DB)7') 

C 

44 FORMAT(1HOf 'DO YOU WANT PRINTOUT TO THE TERMINAL OR TO THE ' 
♦f'LINE PRINTER 'f 

♦/IXf '(TERM OR LPRNT)?') 

C 

45 FORMAT (IHOf 'SPECIFY THE TIME TAG OF EACH FRAME WHICH' f 
♦' IS TO BE PRINTED OUT.'f 

♦/IXf 'INDICATE YOU ARE FINISHED BY ENTERING A -1.0' f 
♦/IXf '(INCLUDE DECIMAL POINT WHEN ENTERING NUMBERS).') 

C 

46 FORMAT (IHOf 14 f' RECORDS HAVE BEEN READ FROM INPUT TAPE') 

C 

47 FORMAT (A8) 

C 

48 FORMAT (IHOf 'DO YOU WANT RAW SYMBOLS TO BE SYNCHRONIZED (1024) 'f 
♦'IN A DATA FRAME OR 'f 

♦/IXf 'THOSE (1025) CORRESPONDING TO A TAPE RECORD (SYNC OR TREC)?' 
C 

49 FORMAT( IHOf 5Xf 'FRAMES DECODED USING QUICK-LOOK ALGORITHM.') 

C 

51 FORMAT (1H0f5Xf 'FRAMES DECODED USING FANO ALGORITHM.') 

C 

52 FORMAT (IHOfSXf 'SYNC WORD LOCATION=' f 15) 

C 

53 F0RMAT(1H0f5Xf 'SIGNAL-TO-NOISE RATIO ( DB ) = ' f G17 . 10 ) 

C 




I 


lei 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

210 

211 

212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 
24F: 


54 F0RMAT<lH0f5X/' PLAYBACK DAY OF YEAR- 'fI3r/) 

C 

55 F0RMAT(lH0»5Xf 'FRAME N0« ' f 15»2Xf 'SYNC 'r 

^ * 'WORD LOCATION*' r 15) 

C SET DEVICE NUMBERS FOR FORTRAN I/O, OF E’OOR QUAUtY 

C IINJ INPUT FROM TERMINAL, 

C IMES: OUTPUT TO TERMINAL, 

C IPRINTtPRINT OUTPUT, 

C ITAPEi: INPUT FROM TAPE, 

C ITAPEOt OUTPUT TO TAPE, 

C IDISK: INPUT AND OUTPUT TO DIRECT ACCESS DEVICE, 

C 

IPASS-0 
50 IIN-1 
IMES-2 
IPRINT«3 
ITAPEI-4 
ITAPEO-5 
IDISK»6 
C 

C FIND OUT WHAT OPTIONS USER WANTS AND SET THE 
C APPROPRIATE FLAGS, 

C 

C FIND OUT WHAT DIRECTION DIGITAL DATA ON TAPE IS STORED IN, 

C 

IFdPASS ,GT, 0) GO TO 200 
100 WRITEdMESf 10) 

X Xrtf / jLJLixr^ 

IF<IDIR,NE,IFWDB,AND,1DIR.NE.IREVB) GO TO 100 
C 

C FIND OUT WHETHER PRINT OR TAPE OUTPUT IS DESIRED, 

C 

200 WRITEdMESf 13) 

READdINf20) lOUT 

IFdOUT,NE,ITP< AND,IOUT,NE,IPR) GO TO 200 
IOPT=IDECA 
220 WRITE <IMESf 44) 

READ d IN f 47) OUTOPT 

IF<0UT0PT,NE,01 ,AND. OUTOPT, NE. 02) GO TO 220 
IF < OUTOPT, EQ, 01) IPRINT=IMES 
C 

C FIND OUT ABOUT PHASE, 

C 

250 WRITEdMESf 43) 

READdlNf20) IPHASE 

IFdPHASE,NE,IDDD,AND,IPHASE,NE,IDBAR) GO TO 250 
IFdOUT,EQ,ITP) GO TO 400 
C 

C SINCE PRINTOUT WAS CHOSEN f FIND OUT IF USER 
C WANTS DATA DECODED OR JUST WANTS TO SEE RAW SYMBOLS, 

C 

300 WRITEdMESf 11) 

READdINf20) lOPT 

IFdOPT,NE.IDECA.AND,IOPT.NE,IRAWA) GO TO 300 
IFdOPT.EQ.IRAUA) GO TO 350 
GO TO 400 
350 WRITEdMESf 48) 

READ<IINf20> IRWD 

IFdRWD.NE,IRTF .AND. IRWD . NE . IRDF ) GO TO 350 
GO TO 600 
C 

C THE DATA IS TO BE DECODED. FIND OUT WHETHER 
C TO USE FANO OR QUICK-LOOK DECODER. 

C 


W'^’. • 

pfl'^^R CjUAS 87V 

400 URITE(IMESrl2> 

READ<IlNr20) ZDEC 

IF<1DEC*NE.IFANA*AND«IDEC»NE.ICILAA) GO TO 400 
IF(IDEC,EQ*IOLAA) GO TO 500 
450 WRITE <IMESf 38) 

REA0<IlN»39rERR*45O) ICOMPT 
IFCICOMPT *EQ* 0) GO TO 460 

IF<IC0MPT.LT«512 .0R» ICOMPT »GT, 10000000) GO TO 450 
GO TO 500 
460 ICOMPT-100000 

SEE IF ALL FRAMES ARE TO BE PROCESSED f OR IF ONLY A 
SELECTED NUMBER OF FRAMES ARE DESIRED* 

500 CONTINUE 
600 WRITE(IMES»15> 

READ < I IN r 20) IANS 

IF(IANS.NE.IYES*AND.IANS*NE*INO) GO TO 600 

IF ALL THE FRAMES ARE TO BE PRINTED OUTf 
SET AN INDICATOR AND MOVE ON* 

IFRAME(1)»-1 

IF<IANS*EQ*IYES) GO TO 1650 

ONLY SELECTED FRAMES ARE TO BE PRINTED OUT. 

FIND OUT IF USER WANTS TO SELECT FRAME BY 
TIME TAG OR FRAME NUMBER* 

700 WRITE(IMESfl6) 

READ < I IN » 20) ITYPE 

IF<ITYPE*NE*ITMA*AND*ITYPE.NE*INMBA) GO TO 700 
IF<ITYPE*EQ*ITMA) GO TO 1100 

USER WANTS TO SELECT FRAMES BY NUMBER. 

: READ IN THE DESIRED FRAME NUMBERS. 

WRITEdMESf 14) 

DO 900 1=1 f 99 

800 READ<IINr22rERR=a00) IFRAME(I) 

IFdFRAMEd) .EQ*0) IFRAMEd)=-l 
IFdFRAMEd) .EQ.-l) GO TO 1000 
900 CONTINUE 

MAKE SURE THERE'S A -1 INDICATING LAST FRAME. 

THEN MOVE ON* 

1000 IFRAME<100)=-1 
GO TO 1650 

USER WANTS TO SPECIFY FRAME BY TIME. ENTER THE TIMES. 

1100 WRITEdMESf45) 

DO 1500 1=1 »99 

: GET HOURS. 


1200 WRITE<IMES»23) 

READ d I N » 26 » ERR= 1200 ) SFR AME < I f 1 ) 

IF<SFRAME<I>1) .EQ.-l. 0) GO TO 1600 

IF<SFRAMEdf 1) .LT.O.O.OR.SFRAMEdd) .GT.24.0)G0T0 1200 




iitsi) 


307 

C 

308 

C GET 1 

309 

C 

310 

311 

312 

313 

1300 1 

1 

314 

C 

315 

C GET 

316 

C 

317 

318 

319 

320 

1400 

321 

1500 

322 

C 

323 

C MAKE 

324 

C THEN 

325 

C 

326 

1600 

327 

C 

328 

C GET 

329 

C 

330 

1650 

331 

332 

333 

1700 

334 

C 

335 

C GET 

336 

C 

337 

338 

339 

340 

341 

342 

343 

1800 

344 

C 

345 

C GET 

346 

C 

347 

348 

1900 

349 

350 

C 

351 

c 

352 

C GET 

353 

354 

355 

C 

356 

C 

357 

C ALL 

358 

C STAR 

359 

C 

360 

C FIRS 

361 

C 

362 

363 

364 

365 

366 

367 

368 

369 

1950 


^“ni, 

OF pocn-i QiiW.UY 


READ < I IN y 26 r ERR-1300 ) SFRAME < I r 2 ) 
IF<SFRAME<Iy2).EQ.-1.0) 00 TO 1600 

IF<SFRAME<Iy2) .LT* 0 . 0 .OR .SFRAMEC y 2) ,GT,60.0)00T0 1300 


READ < I IN y 26 y ERR-1 400 ) SFRAME < I y 3 > 

IF<SFRAME<Iy3) .EQ.-l.O) GO TO 1600 

IF(SFRAME<Iy3) .LT.O.O.OR.SFRAME( I r3) .GT. 60.0) GOTO 1400 


READ<IINy40yERR=1700) IDSS 
IF<IDSS.NE.14.AND.IDSS.NE.43) GO TO 1700 


READtIINy27yERR=1800) lASCyIND 
IFdASC.NE.ISPA.AND.IASC.NE.ILPA) GO TO 1800 
IPBID=11 

IF<IASC.EQ.ILPA) GO TO 1900 
IF<IND.LT.1.0R.IND.GT.3) GO TO 1800 
IPBID=7+IND 


URITEdMESyl?) 
READ<IINy21) ISEQ 


WRITE (IMESy 33) 
READ<IINy21) lYEAR 


IF(IPASS.GT.O) GO TO 2400 
CALL DEFINE<6y5000yl042yIV) 
IF(IV.NE.l) GO TO 5200 
DO 2000 I=ly5000 
READ<ITAPEIy2SyEND=2100) IFR 
IV=I 

URITEdDISK) IFR 


370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 

393 

394 

395 

396 

397 

398 

399 

400 

401 

402 

403 

404 

405 

406 

407 

408 

409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 


C 


OF ' ■ ' ’ ‘ 


2000 CONTINUE 
C 

C SET NUMBER OF RECORDS TO MAX. SHOULD NEVER SET HERE. 

C 

NREC-5000 
GO TO 2200 
C 

C SET THE NUMBER OF RECORDS. 

C 

2100 NREC»I-1 
C 

C SEE IF THE DATA IS RECORDED FORWARDS OR BACKWARDS. 

C 

2200 WRITE <IMESr 46) NREC 

IFdDIR.EQ.IFWDB) GO TO 2400 
C 

C DATA ON TAPE IS IN REVERSE TIME ORDER. REVERSE THE TAPE 
C RECORDS AND SYMBOLS IN EACH RECORD. 

C 

NREC2-<NREC+l)/2 
DO 2300 Ial»NREC2 
NVAL»NREC-I+1 
IV-I 

READ<IDISK) IFR 
IV-NVAL 

READdDISK) IFRl 
CALL REVSlM<IFRj 
CALL REVSIM<IFR1) 

IV=I 

WRITE(IDISK) IFRl 
IV»NVAL 

WRITE(IDISK) IFR 
2300 CONTINUE 
C 

C PRINT THE HEADER FOR THE OUTPUT LISTING. 

C 

2400 WRITE<IPRINTf29) 

IGDIP=IPBID“7 

GO TO <2500r260072700f2800) f IGDIP 
2500 WRITE<IPRINTf30> ISPIA 
IVAL6=I8A 
GO TO 2900 

2600 URITE<IPRINT»30) ISP2A 
IVAL6=I9A 
GO TO 2900 

2700 WRITE<IPRINTf30) ISP3A 
IVAL6=I10A 
GO TO 2900 

2800 WRITE <IPRIMT» 30) ILPIA 
IVAL6=I11A 
C 

C PRINTOUT DSS NUMBER. 

C 

2900 DSS=DSS14 

IF<IDSS.Ea,43) DSS=DSS43 
WRITE<IPRINT.31) DSS 
C 

C PRINTOUT DATA START TIME. THE DATA START TIME 
C WILL BE CONTAINED IN THE FIRST DATA RECORD. SO 
C READ IN THE FIRST RECORD AND GET THE DATA. 

C 


L 


431 

432 

433 

434 

435 

436 

437 

438 
43? 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 

451 

452 

453 

454 

455 

456 

457 

458 

459 

460 

461 

462 

463 

464 

465 

466 

467 

468 

469 

470 

471 

472 

473 

474 

475 

476 

477 

478 

479 

480 

481 

482 

483 

484 

485 

486 

487 

488 

489 

490 


c 


m 

PAGE S3 
OF POOR QUAUI'l Y 

IV*1 

READdDISK) IFR 
C 

C GET PLAYBACK DAY OF YEAR. 

C 

ISDAY«IFR<3)/2 
ISEC-HAND< IVALf 131071 > 

IHOUR-ISEC/3600 
I SEC»I SEC- < IH0UR:tc3600 ) 

IMIN-ISEC/60 

ISEC»ISEC-<IMlN!lc60) 

IHIL»IFR(5) 

SEC0ND-FL0AT< ISEO + . 001#FLOAT< IMIL) 

WRITE<IPRINTf54> ISDAY 
WRITE ( IPRINT r 32 ) IHOUR f IMIN f SECOND 
C 

C SEE IF PRINT OUTPUT IS REQUESTED. IF IT ISN'Tf GO 
r TO TAPE PROCESSING PORTION OF PROGRAM IF ALL FRAMES 
C ARE TO BE PROCESSED. 

C 

II-l 

NFRAMS»0 

NDELT»0 

NSYMBT-0 

NSCT-0 

IF<IQUt.ES. ITP .ANu. lANS.EQ. lYfcS) GO TO 3700 
C 

C FIND OUT IF ALL FRAMES ARE TO BE PROCESSED. 

C 

IFdANS.EQ.IYES) GO TO 3950 
C 

C ONLY SELECTED FRAMES ARE TO BE PROCESSED. PROCESS EACH FRAME. 
C 

C IF TAPE OUTPUT OF SELECTED FRAMES IS DESIRED f FIRST WRITE 
C OUT FILE HEADER. 

C 

IF(IOUT.EQ.ITP) GO TO 3700 
2950 K=1 

IFCIOUT.EQ.ITP .AND. IDEC.EQ. IQLAA) WRITE( IPRINT. 49) 
IFdOUT.EQ.ITP .AND. IDEC.NE. IQLAA) WRITE< IPRINT. 51 ) 

C 

C FIRST SEE HOW FRAME WAS SELECTED. 

C 

3000 IF<ITYPE.EQ.INMBA) GO TO 3100 
C 

C USER WANTS TO FIND FRAME BY TIME REFERENCE. 

C 

IF<SFRAME<II.i) .EQ.-1.0.0R.SFRAME(II.2) .EQ.-l.O.OR. 
;|CSFRAME<II.3) .EQ.-l.O) GO TO 3600 
C 

C CALL A ROUTINE TO FIND THE FRAME NUMBER CORRESPONDING TO THE 
C TIME TAG SPECIFIED, 

C 

CALL FRFINDdl.SFRAME.IDISK.IFR.NREC.NFR.IU) 

IF<NFR.LT.l) GO TO 3400 
NFRAMS=NFRAMS+1 
GO TO 3200 
C 

C USER SELECTED FRAME BY NUMBER. 

C 


3100 NFR*IFRAME<II> 

IF<NFR.EQ»-1) 00 TO 3600 
IF<NFR*LT.1.0R,NFR.0T*NREC> 00 TO 3400 
NFRAMS-NFRAMS+1 


aicr'i PAGE « 
rtF PC-ON QUAUfV 


C NOW HAVE THE FRAME NUMBER. SEE IF RAW SYMBOL OUTPUT IS DESIRED. 

C « 

3200 IF(IOPT.EO.IRAUA) GO TO 3500 
C » 

C USSR WANTS DATA DECODED. FIRST FIND THE SYNC WORD AND GET THE FRAME. 

C 

CALL SYNC ( NFR r NREC » "DISK » IFRB 1 1 WORK r IRETA r IQLA f 
*lpESNOrIFFRTpISYCWfIV) 

IF<IFFRT.QT.O) 00 TO 3250 

NDELT-NDELT+1 

GO TO 3400 ' 

C 

C SEE WHAT KIND OF DECODER USER WANTS. 

C 

3250 IDFFLQ<1)»-1 

IF(IOUT.EQ.ITP) WRITE (IPRINTf 55) NFRrISYCW 
. IF(IOUT.EQ.ITP) WRITE<IPRINTf53) ESNOCl) 

. IFdDEC.EO.IQLAA) GO TO 3300 
C 

C USER WANTS FANO ALGORITHM USEDr SO CALL ROUTINF TO DO IT. 

C 

ITCT<1>=IC0MPT 

CALL FAN0<ESN0<1) fIRETApIQLArITCT<l) f ACB(l) pNSC(I) f IDFFLG<1 ) pIOUT 
NDELT=NDELT+1 

. IF(IDFFLG<1) .NE.O) GO TO 3300 
''*'NDELT=sNDELT-1 

NSYMBT-NSYMBT+1024 

NSCTaNSCT+NSC<l> 


C DATA IS DECODED PRINT IT OUT. 


526 . C 

) 527 3300 IF(IOUT.EQ.ITP) GO TO 6000 

j 528 : CALL FRPRNT<NFRfIQLArACBfNSCpIDFFLQf IDISKf IFRp 

529 - *IDECf If IPRINTf ISYCWfESNOpIV) 

, 530 c; 

531 D GO ON AND GET NEXT SELECTED FRAME. 

532 C* 

533 3400 11=11+1 

534 . GO TO 3000 

535 C 

536 C- RAW SYMBOL OUTPUT REQUESTED. CALL ROUTINE TO DO IT, 

537 C 

538 C FIRST GET S/N RATIO AND SYNCHRONIZE THE (MAPPED) SYMBOLS 

539 C. 

540 3500 CALL SYNC(NFRfNRECf IDISKf IFRBf IWORKf IRETAf IQLAf 

541 #lfESNOf IFFRTf ISYCWy Z V 

542 C* ^ 

543 C. TRANSFER (MAPPED) SYMBOLS AND S/N RATIO TO IWORK ARRAY. 

544 6 FIRST UNCOMPLEMENT THE SECOND SYMBOL OF EACH CHANNEL PAIR, 

545 Cf 

546 ' . DO 3525 1=1 f 1023. 2 

547 J=I+1 

548 IWORKdf l) = IRETA(If 1) 

549 • IWORK(Jf l)=9~IRETA(Jf 1) 

550 3525 CONTINUE 

551 • IW0RK(1025f 1)=ISYCW 

■ - ----^1-.. . « ^ o 1 


552 

553 

554 

555 

556 

557 

558 

559 

560 

561 

562 

563 

564 

565 

566 

567 

568 

569 

570 

571 

572 

573 

574 

575 

576 

577 

578 

579 

580 

581 

582 

583 

584 

585 

586 

587 

588 

589 

590 

591 

592 

593 

594 

595 

596 

597 

598 

599 

600 
601 
602 

603 

604 

605 

606 

607 

608 

609 

610 
611 
612 

613 

614 


C 

C 

C 


CALL RSPRNT < NFR » IDISK f IFR f IPRINT f I WORK r NREC t ESN0 1 
JKIRWrirlPHASEf IV) 

GO TO 3400 


ALL DECODING DONE, PRINT SUMMARY BLOCK, 

3600 IF<NFRAMS ,EQ, 0) GO TO 3650 
II»II-1 

•IF<IOUT,EQ,ITP) GO TO 6002 
3650 GO TO 4700 
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C 

C 

C 


TAPE OUTPUT IS REQUESTED, GENERATE AND WRITE OUT FILE HEADER- 

3700 NFRAMS=0 
NDELT=0 
NSYMBT=0 
NSCT=0 
IFLBR<1)=0 
IFLBR<2)=84 
IFLBR<3)=ILMA 
IFLBR<4)=IIDA 
IFLBR<5)=20992 
IFLBR<6)e25 

IF<IDSS,EQ,43) IFLBR(6)=134 

IVAL5=IUAL6 

IFLBR<9)=0 

IFLBR < 10 ) =MLSR2 < ISEQ r 8 ) 

IFLBR<10)=IFLBR<10)~192 

IFULL=ISEQ 

IFLBR < 1 1 ) = < HAND < IFULL r 255 ) -1 92 ) *256 
IF=70 

IF(IDIR,EQ,IREVB) IF=82 

IFLBR(11)=IFLBR<11)+IF 

IFLBR<12)=IWD 

IF<inXR,EQ,IREVB) IFLBR< 12)=IEV 
IFLBR (13 ) = < MLSR2 < lYEAR , 8 ) -192 ) +48*256 
IFULL=IYEAR 

IFLBR < 14 )= (MAND< IFULL r 255 >-192 >*256 
DO 3800 III=15»23 
IFLBR<III)=0 
3800 CONTINUE 

GET DATA START TIMEfDAY OF YEAR, 

IV=1 

READ(IDISK) IFR 
ISEC=MAND( IVALf 131071 ) 

IYEARS=IFR<3)/2 
IFLBR ( 24 )=I YEARS 
IMIL=IFR<5) 

IHSEC=IMIL/10+ISEC*100 
GET DATA STOP TIME, 

IV=NREC 

READCIDISK) IFR 
IYEARS=IFR<3>/2 
IFLBR (25) =IYEARS/256 
IFLBR < 26 ) =MLSL2 < I YEARS , 8 ) 

IFULL=IHSEC 

IFLBR < 26 ) = IHALF < 1 ) +IFLBR < 26 ) 

IFLBR(27)=IHALF(2) 

ISEC=MAND(IVALf 131071) 


C 

C 

c 


c 

c 

c 






615 

616 

617 

618 

619 

620 
621 
622 

623 

624 

625 

626 

627 

628 

629 

630 

631 

632 

633 

634 

635 

636 

637 

638 

639 

640 

641 

642 

643 

644 

645 

646 

647 

648 

649 

650 

651 

652 

653 

654 

655 

656 

657 

658 

659 

660 
661 
662 
663 


IMIL=IFR<5> ' 

IHSEC=iriIL/10+ISEC3K100 . .. . • 

IFULL»IHSEC 
IFULL.=MLSL<IFULLr8) 

IFLBR(28)=IHALF(1) 

IFLBR(29)«IHALF<2> 

DO 3900 IHa30f608 
IFLBR(III)«0 
3900 CONTINUE 

WRITE OUT THE FILE LABEL RECORD, 

WRITE<ITAPE0»37) IFLBR 

SET UP THE COMMON PART OF THE DATA BLOCK. 

IFLBR(1)=1 
IFLBR(2)=1216 
IFLBR(3)=150 
IFLBR<4)=150 
IFLBR<5)=8 
IFLBR<6)=0 
IFLBR<7)=0 
IFLBR (8>=C 


*i PASS ts 


SAUTY 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


IF TAPE OUTPUT OF SELECTED FRAMES IS DESIRED RETURN 
"0 DECODE THOSE FRAMES, 

IF<IAN3,EQ.IN0) GO TO 2950 


C 

c 

c 


IF ALL FRAMES ARE DESIRED START DECODING r 
EIGHT FRAMES AT A TIME, 

3950 IFaOPT,EQ,IRAWA) GO TO 5300 
Jal 

4000 CALL SYNC( JfNRECf IDISKf IFRBf lUORKf IRETAf IQLAf 
♦SrESNOrlFFRTf ISYCWrIV) 

IF(IFFRT,LE,0) GO TO 4400 
NFRAMS=NFRAMS+IFFRT 

IF QUICK-LOOK CODE IS BEING USEDf DON'T CALL FANO DECODER. 

IFdDEC ,NE, IQLAA) GO TO 4030 
IF<IOUT.EQ,ITP) WRITE(IPRINTf49) 

IF<IOUT,EQ,ITP) WRITE(IPRINTf52) ISYCW 
DO 4020 I=ifIFFRT 

IF(IOUT,EQ,ITP) WRITE<IPRINTf53) ESN0<I) 

4020 CONTINUE 

4030 IF<IDEC,EQ, IQLAA) GO TO 4150 


664 C 

665 C CALL FANO DECODER TO DECODE EACH FRAME, 

666 C 

667 IFdOUT.EQ.ITP) URITE( IPRINTf 51 ) 

668 IF<IOUT,EQ,ITP) WRITE< IPRINT»52> ISYCW 

669 DO 4100 I=lrIFFRT 

670 IF<IOUT.EQ,ITP) WRITE <IPRINTf53) ESNO(I) 

671 ITCT(I)=ICOMPT 

672 NDELT=NDELT+1 

673 CALL FANO<ESNO(I)fIRETA<lfI) rlQLAdfl) f ITCT(I) fACB<I) r 

674 *NSC(I) f IDFFLG<I)rIOUT) 

675 IFdDFFLGd) .NE.O) GO TO 4100 

676 NDELT=NDELT-1 


• t 


677 

678 

679 

680 
681 
682 

683 

684 

685 

686 

687 

688 

689 

690 

691 

692 

693 

694 

695 

696 

697 

698 

699 

700 

701 

702 

703 

704 
70d 

706 

707 

708 

709 

710 
711. 

712 

713 

714 

715 

716 

717 

718 

719 

720 

721 

722 

723 

724 

725 

726 

727 

728 

729 

730 

731 

732 

733 

734 

735 

736 

737 

738 

739 


NSYMBT-NSYMBT+1024 
NSCT*NSCT+NSC<I) 
4109 CONTINUE 


O.;' W®. '3 

OF F 




C SEE IF TAPE OUTPUT IS DESIRED* 

C 

4150 IF<I0UT*EQ*1TP) 60 TO 4175 
C 

C PRINTOUT IS DESIRED. CALL ROUTINE TO DO THE PRINTING. 
C 


CALL FRPRNT ( J f I QLA r ACB r NSC r IDFFL6 r I D I SK » I FR » 
5KIDEC f IFFRT r IPRINT r ISYCW r ESNO r IV ) 

60 TO 4400 


C 

C TAPE OUTPUT REQUESTED. WRITE OUT THE DATA BLOCKS. 
C 

4175 DO 4200 I«lr IFFRT 


INDEX«<I-1)>IC75 

IARG=9+INDEX 

IFLBR<IAR6)«25206 

IARG=10+INDEX 

IFLBR<IARG)=25 

IF<IDSS.EQ.43> IFLBR< IARG)=134 
IVALLL=IFLBR(IARG> 

IFLBR ( I ARC ) =MOR ( I VALLL r 9984 ) 
IFULL=9728+IPBID 

IFdPBID .LT. 11) IFULL=8960+IPBID 


IFULL= I FULL+ 1 0665984 
IFULLsMLSL < IFULL 7 8 ) 


IARG=12+INDEX 
IFLBR tlARG ) alHALF ( 1 ) 

IARG=13+INDEX 

IFLBR(IARG)=IHALF(2) 

IARG=11+INDEX 

IFLBR ( I ARG ) =MOR ( MLSL ( 1 70 f 8 > f 222 ) 

IV=J+I 

READ<IDISK) IFR 
ISEC=MAND< IVALi 131071 ) 

IHIL=IFR<5) 

IHSEC=IMIL/10+ISEC*100 

IFULL=IHSEC 

IARG=13+INDEX 

IFLBR < I ARG ) = IFLBR < lARG ) + IHALF < 1 > 

IARG«14+INDEX 

IFLBR<IARG)=IHALF(2) 

IYEARS=IFR(3)/2 
JDAYH=IYEARS/100 
' IDAYT= < IYEARS-IDAYH5K100 ) / 10 
IDAYU=IYEARS-IDAYH3lcl00-IDAYT5K10 
IARG=15+INDEX 

IFLBR<IARG)«IDAYH*4096+IDAYT*256+IDAYU*16 

IAR6al6+INDEX 

IFLBR< IARG)«IMIL/100+I*256 
IDAYH=128 

IFdDFFLGd) .NE.ON IDAYH=130 
IAR6=17+INDEX 

IFLBR < I ARG > =MLSL2 ( IDAYH r 8 ) 


IAR3=18+INDEX 
IFLBR d ARG )=0 
IARG=19+INDEX 
IFLBRdARG) = IFR(8) 
IARG=20+INDEX 
IFLBRdARG)fIFR<7) 


740 

741 

742 

743 

744 

745 

746 

747 

748 

749 

750 

751 

752 

753 

754 

755 

756 

757 

758 

759 

760 

761 

762 

763 

764 

765 

766 

767 

768 

769 

770 

771 

772 

773 

774 

775 

776 

777 

778 

779 

780 

781 

782 
• 783 

784 

785 

786 

787 

788 

789 

790 

791 

792 

793 

794 

795 

796 

797 

798 

799 

800 
801 


IARG=*21+INDEX 

CALL DPACK < ICJLA < 1 f I > » IFLBR < I ARG ) ) 

IARG=53+INDEX 

IFLBR<IARG)=NSC(I) 

. IARG=54+INDEX 
IFLBR ( I ARG ) =ITCT ( I ) /64 
DO 4200 I I 1*55 r 83 
lARGalll+INDEX 
IFLBR<IARG)=0 
4200 CONTINUE 
C 

C SEE IF FILLER BLOCKS HAVE TO BE INSERTED* 
C 

IF<IFFRT.EQ.8) GO TO 4350 
C 

C WRITE UP FILLER BLOCKS. 

C 

IFRT1=IFFRT+1 

DO 4300 III=IFRT1»8 

INDEX=<III-1)*75 

IARG*9+INDEX 

IFLBR<IARG)=25206 

IARG=10+INDEX 

IFLBR<IARG)=25 

IFdDSS .EQ. 43) IFLBR< IARG) = 134 
IVALLL=IFLBR<IARG) 

IFLBR < I ARG ) =MOR < I VALLL r 9984 ) 
IARG*11+INDEX 

IFLBR < I ARG ) =MOR < MLSL < 1 70 r 8 ) r 1 90 ) 

DO 4300 1111=1^72 
IARG=11+INDEX+IIII 
IFLBR(IARG)=IFILL 
4300 CONTINUE 

4350 WRITE (ITAPEOr 37) IFLBR 
IFLBR(1)=IFLBR<1)+1 
C 

C PROCESS NEXT GROUP OF FRAMES. 

C 

4400 J=J+8 

IF(J.GT.NREC) GO TO 4500 
GO TO 4000 


OF QUALITY 


c 

C DATA DECODED. NOW WRITE TAPE ENDING IF TAPE OUTPUT REQUESTED. 

C 

4500 IF(IOUT.NE.ITP) GO TO 4700 
IFLBR(1)=IKEX1 
DO 4600 I=2r608 
IFLBR(I)=0 
4600 CONTINUE 

WRITE (ITAPEOf 37) IFLBR 
END FILE ITAPEO 
END FILE ITAPEO 

C 

C WRITE SUMMARY BLOCK FOR THIS PASS THROUGH THE PROGRAM. 


4700 IF(IOPT.EQ.IDECA.AND.IDEC.EQ.IFANA) GO TO 4800 
WRITE(IPRINTf35) NFRAMS 
GO TO 4900 
4800 DELRT=0.0 
SERT=0.0 

IF < NFRAMS . NE . 0 ) DELRT=FLOAT ( NDELT ) /FLOAT < NFRAMS ) 
IF < NSYMBT . NE . 0 ) SERT*FLOAT ( NSCT ) /FLOAT ( NSYMBT ) 


0 


802 

803 

804 

805 

806 

807 

808 

809 

810 
811 
812 

813 

814 

815 

816 

817 

818 

819 

820 
821 
822 

823 

824 

825 

826 

827 

828 

829 

830 

831 

832 

833 

834 

835 

836 

837 

838 

839 

840 

841 

842 

843 

844 

845 

846 

847 

848 

849 

850 

851 

852 

853 

854 

855 

856 

857 

858 

859 

860 
861 
862 

863 

864 


WRITE ( IPRINT 1 34 ) NFRAMS t NDELT y DELRT r SERT 
4900 WRITE <IMESf 36) 

READ ( I IN r 20) IANS 

IF(IANS«NE.IYES,AND*IANS»NE»INO) GO TO 4900 
IF<IANS«EQ,IYES) GO TO 5000 
IPASS=IPASS+1 
GO TO 50 
5000 STOP 
C 

C ERROR IN ACCESSING TEMPORARY FILE. 

C 

5200 WRITE<IMESf42) 

STOP 
C 

5300 NFRAMS=0 

DO 5400 I=1»NREC 
NFRAMS=NFRAMS+1 
CALL RSPRNT < I » IDISK r IFR f IPRINT f I WORK » NREC f ESNO r 
!|cIRWD»IPHASE»IV) 

5400 CONTINUE 

GO TO 4700 
C 

C OUTPUT TO TAPE THE FRAMES SELECTED AFTER PROCESSING. 
C 

6000 INDEX=<K-1)*75 
IARG=9+INDEX 
IFLBR(IARG) =25206 
IARG=10+INDEX 
IFLBR(IARG)=25 

IF<IDSS .EQ. 43) IFLBR ( lARG ) =134 
IVALLL=IFLBR<IARG) 

IFLBR ( I ARG ) =MOR < I VALLL f 9984 ) 

IFULL=9728+IPBID 

IFdPBID .LT. 11) IFULL=8960+IPBID 
IFUl. L=IFULL+10665984 
IFULL=MLSL(IFULLf8) 

IARG=12+INDEX 
IFLBR <1 ARG ) =IHALF ( 1 ) 

IARG=13+INDEX 
IFLBR < I ARG ) =IHALF ( 2 ) 

IARG=11+INDEX 

IFLBR < I ARG ) =MOR < MLSL ( 170 f 8 ) r 222 ) 

IV=NFR+1 
READdDISK) IFR 
ISEC=MAND<IOALr 131071) 

IMIL=IFR(5) 

IHSEC=IMIL/10+ISEC*100 

IFULL=IHSEC 

IARG=13+INDEX 

IFLBR ( I ARG ) =IFLBR < lARG ) +IHALF ( 1 ) 

IARG=14+INDEX 
IFLBR < I ARG ) =IHALF ( 2 ) 

IYEARS=IFR(3)/2 

IDAYH=IYEARS/100 

IDAYT=(IYEARS-IDAYH)K100)/10 

IDAYU=IYEARS-IDAYH>K100~IDAYT:«?10 

IARG=15+INDEX 

IFLBR(IARG)=IDAYH5lc4096+IDAYT*256+iriAYUHa6 

IARG=16+INDEX 

IFLBR<IARG)=IMIL/100+K*256 

IDAYH=128 

IFdDFFLGd) .NE.O) IDAYH=130 
IARG=17+INDEX 


OH' 

OF 




XITY 


865 

866 

867 

868 

869 

870 

871 

872 

873 

874 

875 

876 

877 

878 

879 

880 
881 
882 

883 

884 

885 

886 

887 

888 
889 
390 

891 

892 

893 

894 

895 

896 

897 

898 

899 

900 

901 

902 

903 

904 

905 

906 

907 

908 

909 

910 

911 

912 

913 

914 

915 

916 

917 

918 

919 

920 

921 

922 

923 

924 

925 


yi- POUR QOALSTY 


c 

c 

c 


IFLBR < I ARC ) =»MLSL2 < IDAYH r 8 ) 

IARG»18+INDEX 
IFLBR<IARG)«0 
IARG»19+INDEX 
■ IFLBR<IARG)aIFR<8) 

IARG»20+INDEX 

IFLBR<IARG)»IFR(7) 

IARG«21+INDEX 

CALL BRACK ( IQLA < 1 r 1 ) f IFLBR < I ARG ) > 

IARG*53+INDEX 

1FLBR<IARG)=NSC(1) 

IARG»54+INDEX 
IFLBR < I ARG ) => ITCT < 1) /64 
DO 6001 III»55r83 
lARQalll+INDEX 
IFLBR (I ARG )=0 
6001 CONTINUE 

CHECK IF DONE AND SEE IF FILLER BLOCKS HAVE TO BE INSERTED* 


C 

C 

C 


K=K+1 

IF<K *EQ* 9) GO TO 6004 
KK=II+1 

IF<SFRAME(KK»1) *NE*-1.0 *AND* SFRAME(KKf 2) *NE*-1 .0 
#SFRAME<KKf3)*NE.-1.0 .AND* IFRAME<KK) *NE*-1 ) 60 TO 

WRITE UP FILLER BLOCKS* 


.AND* 

3400 


6002 DO 6003 III=K»8 
INDEX=<III-1)*75 
IARG=9+INDEX 
IFLBR < I ARG >=25206 
IARG=10+INDEX 
IFLBR<IARG)=25 

IFdDSS *EQ* 43) IFLBR< IARG) = 134 
IVALLL=IFLBR<IARG) 

IFLBR < lARG ) =MOR ( I VALLL r 9984 ) 

IARG=11+INDEX 

IFLBR < I ARG ) =MOR < MLSL ( 170 » 8 ) f 190 ) 

DO 6003 IIII=1»72 

IARG=11+INDEX+IIII 

IFLBR<IARG)=IFILL 

6003 CONTINUE 

6004 WRITE (ITAPEOf 37) IFLBR 
IFLBR<1)=IFLBR(1)+1 
K=1 

KK=II+1 

IF<SFRAME<KKrl> .NE.-l.O .AND. SFRAME(KKf2) .NE.-l *0 .AND* 
:»cSFRAME(KKf3) .NE.-l.O .AND. IFRAME(KK) .NE.-l ) GO TO 3400 


C 

C 

C 


IF ALL SELECTED FRAMES ARE PROCESSED WRITE TAPE ENDING. 


IFLBR<1)=IHEX1 
DO 6005 I=2f608 
IFLBR<I)=0 
6005 CONTINUE 

WRITE<ITAPE0f37) 
END FILE ITAPEO 
END FILE ITAPEO 
GO TO 4700 
END 


IFLBR 




Q. 


G 


1 

2 

3 

4 
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SUBROUTINE REMSIH<IFR) 

C 

C 

C THIS SUBROUTINE REVERSES THE SYMBOLS IN AN INPUT BATA 
C RECORD. 

C 

C DIMENSION THE ARRAY AND DECLARE THE HALFWORD INTEGERS. 

C 

INTEGER!»C2 IFR ( 521 ) f I VALl f I VAL2 r IFLIP 
C 

C REVERSE EACH WORD AND REVERSE THE SYMBOLS IN EACH WORD. 
C 

DO 100 I»lr256 

IARG=521-I+1 

IVALl=IFLIP<IFR(I+9)> 

IVAL2=IFLIP<IFR<IARG) ) 

IFR(I+9)*IVAL2 
IFR<IARG)«=IVAL1 
100 CONTINUE 
RETURN 
END 

FUNCTION IHARD(Il) 

C 

C 

C THIS FUNCTION BENEKATES A HARD DECISION. 

C 

INTEGER*2 II 
IHARD=1 

IF<I1.LE.4) IHARD=0 

RETURN 

END 

INTEGER FUNCTION MAP:«c2(Il) 

COMMON IPHASE 
DATA IDB/'DB V 
C 
C 

C FUNCTION TO MAP SOFT-DECISIONS. 

C 

INTEGER:|C2 II 

IFdPHASE.EQ.IDB) GO TO 100 
C 

C 'DATA' PHASE REFERENCE. 

C 

MAP=12-I1 

IF<I1.LT.4) MAP=4-I1 
RETURN 
C 
C 

C 'DATA-BAR' PHASE REFERENCE. 

C 

100 MAP=5+I1 

IF(I1.GT.3) MAP=Il-3 

RETURN 

END 

SUBROUTINE SYNC < JS r NREC r IDISK » IFRB r I WORK r IRETA r 
JtclOLA r NFRMS r ESNO r IFFRT r ISYCW r I V ) 

C 

C 

C THIS SUBROUTINE LOCATES THE SYNC WORD. 


60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 


C 

C 

C 


C 

C 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


DECLARE VARIABLES AND DIMENSION ARRAYS. 

INTEQERHC2 IFRB<521 r 9> f IW0RK<1025f 9) r IRETA< 1024r 8) r 
*IQLA<512r9)f ISU0RD<23) 

Dimension ESNO<s>fipos<2) risYNc< 2 ) 

DATA ISWORD/lflflflrlfOfOfO/lflfOfOfOflf 
*OrlrOflrOfOflfOfO/ 

SEE WHAT WE HAVE TO WORK WITH. 

IFFRT=NFRMS 

IF < NFRMS+ JS . GT . NREC ) IFFRT=NREC- JS 
IFCIFFRT.LE.O) GO TO 1100 
IFRT1=IFFRT+1 
J JS 


READ THE FRAMES IN. 

DO 100 I=lfIFRTl 

INDEX=JJ+I-1 

IV=INDEX 

READ<IDISK) aFRB(J»I> f J=lf521) 

GET THE SIGNAL-TO-NOISE RATIO. 

IESN0=IFRB(8?I) 

ESNO < I ) =FLOAT ( lESNO > /128 . 0 
100 CONTINUE 

EXPAND THE FRAMES. 

DO 200 I=lrIFRTl 

CALL UNPCKCIFRB(lrl) rlWORK(lrl) ) 
200 CONTINUE 

GO THROUGH BOTH SYMBOL PAIRS. 


U£- 


QIAUIY 


98 

99 
100 
101 
102 
103 


ISYNC(l)=-25 
ISYNC<2)=-25 
DO 500 1=1.2 
C 

C COMBINE THE SYMBOLS. 
C 


104 


DO 300 J=1.IFRT1 


105 

106 
107 


DO 300 K=1.512 

INDEX=2*K+I-2 

INDEX1=INDEX+1 


108 

109 


IVAL=IHARD<IWORK(INDEX. J) ) 
IVAL1=1“IHARD(IW0RK(INDEX1»J) ) 


110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 
121 


IQLA ( K > J ) =MXQR <1 VAL . I VALl ) 
300 CONTINUE 
C 

C NOW SEARCH FOR THE SYNC WORD. 

C 

DO 500 J=1.512 
IC0L=0 

DO 400 K=1.IFRT1 
DO 400 L=lf23 
LL=J+L-1 

IF(LL.GT.512) LL=LL-512 
I0R=1 



122 



IF < ISWORD < L ) ♦ NE ♦ IQLA < LL r K ) ) IOR—1 


123 



ICOL*ICOL+IOR 


124 


400 

CONTINUE 

% 

125 



IF<ICOL.LT.ISYNC<I>) GO TO 500 


126 



IPOS<I)=J 


127 



ISYNC(I) = ICOL 

1 

128 


500 

CONTINUE ‘ • i VV-. 


129 

C 


OF OL ' • 


130 

C 

LOCATE SYNC WORD. 


131 

c 



3 

132 



lal 


133 



IF<ISYNC<2).GT.ISYNC<1>) 1=2 


134 



LL»IP0S<I)+24-l 

9 

135 



IF<LL.GT.S12) LL=LL-512 


136 



ISYCUaLL 


137 



LL=LL^c2-2+I 

• 

1 

138 



LLL=LL 


139 

c 




140 

c 

NOU 

SET UP THE RETURN ARRAY WITH THE SOFT-DECISIONS 


141 

c 



• 

142 



DO 800 I=lfIFFRT 


143 



DO 700 J«lfl024 


144 



JJ=LL+J-1' 

\ 

145 



IF<JJ.NE.1025) GO TO 600 

i 

146 



JJ=JJ+1 

I 

147 



LL=LL+1 


148 


600 

IRcTA< Jf I )=IWORK< JJ» I) 

1 

149 


700 

CONTINUE 

J 

150 



LL=LLL 


151 


800 

CONTINUE 

1 

152 

c 



J 

153 

c 

NOU 

GO BACK AND CALCULATE THE HARD-DECISIONS AND 


154 

c 

THE 

OUICK-LOOK OUTPUTS. 


155 

c 


m 


156 



DO 1000 I=lfIFFRT 


157 



DO 900 J=lr511 


158 



INDEX=2*J+1 


159 



INDEX1=INDEX+1 

, 

160 



I VAL=IHARD < IRETA < INDEX r I ) ) 


161 



IVALl=l-IHARD(IRETA<INDEXlf I) ) 


162 



IRETA< INDEXl r I )=9-IRETA ( INDEXl r I ) 


163 



IQLA < J » I ) =MXOR ( IVAL f IVALl ) 

* 

164 


900 

CONTINUE 


165 



IRETA < 2 f I ) =9-IRETA < 2 » I ) 


166 



IQLA(512»I)=1 

i 

167 

1000 

CONTINUE 


168 

1100 

RETURN 


169 



END 


170 



SUBROUTINE FRPRNT(NFRf IQLAf ACB»NSCf IDFFLGr IDISK 


171 


*IDECfICNTfIPRINTf ISYCUfESNOflV) 

* 

172 

c 




173 

c 



a 

174 

c 

THIS SUBROUTINE PRINTS OUT FRAMES* 


175 

c 




176 

c 

DECLARE VARIABLES AND DIMENSION ARRAYS. 

• 

177 

c 




178 



INTEGERHC2 IQLA<512f 8 ) f IFR < 521 ) f IFRD ( 2 ) 


179 



DIMENSION ACB(8) r NSC< 8 ) r IDFFLG< 8 ) f ESNO ( 8 ) 


180 



EQUIVALENCE ( IVALf IFRD ( 1 ) ) 

a 

181 



DATA IDECA/'FANOV 


182 

c 




183 

c 

FORMAT STATEMENTS USED BY ROUTINE. 


184 

c 




tQ > 


185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 

223 

224 

226 

227 

228 

229 

230 

231 
r>~Kr> 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 


10 F0RMAT<lH0f//6Xf 'PRINTOUT OF DECODED DATA FOR FRAME 
■*'DATA START TIME-'f/lOXf 

5K'H0URS- 'fI2r/10Xf 'MINUTES- '» 

)KI2f/10Xr 'SECONDS- 'fF6.j;> 

11 FORMAT <lH0r/6Xf 'DATA STOP TIME-' r/lOXr 'HOURS- 'rI2r/10Xr 
:K'M1NUTES- ' rI2f/10Xr 'SECONDS- 'rF6*3) 

12 FORMAT (IHOrSXf 'FRAME WAS DECODED SUCCESSFULLYf ' > 

13 FORMAT (IHOfSXf 'FRAME WAS DELETED*') 

14 FORMAT <lH0r5X» 'FRAME WAS DECODED USING QUICK-LOOK 'r 
♦'ALGORITHM*') 

15 FORMAT <1H0»5X» 'FRAME WAS DECODED USING FANO '» 

♦'ALGORITHM*') 

16 FORMAT (lH0f5X»' NUMBER OF SYMBOL ERRORS CORRECTED- '» 

♦ 14) 

17 FORMAT <lHOfSX»' AVERAGE NUMBER OF COMPUTATIONS PER BIT«'rF14*7) 

18 FORMAT UHOfSXf 'DECODED DATA-' r//3Xr64Il »/3Xr64Il f 
♦/3Xf64Ilf/3Xf64Ilf/3X»64Ilf/3Xf64Ilr/3Xf64Il»/3Xf64Il) 

19 FORMATdHOf' TAIL SEQUENCE SHOULD BE- ', 

♦'lllllOvOllOOOlOlOlOOlOOl' ) 

20 F0RMAT<lH0r5Xf' INVALID DECODER PARAMETER*') 

21 FORMAT (IHOfSXr 'SYNC WORD LOCATION-' » IS) 

22 F0RMAT(lH0f5X? 'SIGNAL-TO-NOISE RATIO <DB)-' r G17* 10) 

23 FORMATdH fSXf'C* BITS DECODED/MAX* ♦ COMP* PER FRAME)') 

C 

C LOOP THROUGH EACH FRAME. 

C 

DO 600 I=lfICNT 
C 

C GET DATA START TIME* 

C 

INDEX-NFR+I-1 
IV-INDEX 
READdDISK) IFR 
IFRD<1)=IFR<3) 

IFRD(2)=IFR(4) 

ISEC=MAND( IVALf 131071 ) 

IMIL=IFR(5) 

IHOUR-ISEC/3600 
ISEC=ISEC-<IH0UR^3600) 

IMIN-ISEC/60 
ISEC=ISEC-<IMIN^60) 

SECOND-FLOAT ( ISEC ) + . OOl^FLOAT < IMIL) 

C 

C PRINT OUT DATA START TIME. 

C 

WRITE<IPRINTf 10) INDEX r IHOURr IMINr SECOND 
C 

C GET DATA START TIME. 

C 

INDEXl-INDEX+1 
IV-INDEXl 
READ(IDISK) IFR 
IFRD(1)=IFR(3) 

IFRD<2)=IFR<4) 

ISEC=MAND( IVALf 131071 ) 

IMIL=IFR(5) 

IHOUR-ISEC/3600 
ISEC=ISEC-( IH0UR^3600 ) 

IMIN-ISEC/60 
ISEC=ISEC-<IMIN^60) 

SECOND-FLOATdSEO + .OOl^FLDATdMIL) 


Ij -' 








4 


2A6 

247 

248 


250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 


274 

275 

276 

277 

278 

279 

280 
281 
282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 

302 

303 

304 

305 

306 
307 


C 

C 

C 


C 

C 

C 


C 

C 

C 


C 

C 

C 


C 

C 

C 


c 

c 

c 


PRINT OUT DATA STOP TIME. 

WRITE< IPRINTr 11 > IHOURr IMINr SECOND 
PRINT OUT SYNC WORD LOCATION .,ND S/N RATIO 


WRITE(IPRINT»21) ISYCW 
WRITE<IPRINT»22) ESNO(I) 


CHECK ON DECODER. 

IF<IDEC.EO.IDECA) GO TO 100 
QUICK-LOOK DECODER USED* 


OF PDUfi 


WRITE<IPRINTrl4) 
GO TO 200 


C 

c 

c 


FANO DECODER USED. 

100 WRITE(IPRINT»15) 

SEE IF FRAME WAS DECODED SUCCESSFULLY. 

200 IF<IDEC.EQ.IDECA.AND.IDFFLG<I) .NE.O) GO TO 300 
FRAME WAS DECODED* 


IF(IDEC.NE.IDECA) GO TO 400 
WRITE(IPRINT»12) 

GO TO 400 


C 

C 

C 


FRAME WAS DELETED. 
300 WRITEdPRINTf 13) 


C 

C 

c 

c 


OUTPUT NUMBER OF SYMBOL ERRORS AND AVERAGE NUMBER OF 
COMPUTATIONS PER BIT. 


400 IFdDEC.NE. IDECA) GO TO 500 
WRITEdPRINT»16) NSCd) 
WRITEdPRINTf 17) ACBd) 
IFdDFFLGd) .EQ. 0) GO TO 500 
WRITEdPRINTf 23) 


C 

C 

C 


OUTPUT THE DECODED FRAME. 


500 WRITEdPRINTf 18) (IQLA< Jf I ) f J=1 f 512 ) 

IFdDFFLGd) .EQ.O.OR.IDEC.NE. IDECA) GO TO 600 
IFdDFFLGd) .NE.l) GO TO 550 
WRITEdPRINTf 19) 

GO TO 600 

550 WRITEdPRINTf 20) 

600 CONTINUE 
RETURN 
END 

SUBROUTINE FRFIND < 1 1 f SFRAME f IDISK . IFR f NREC . NFR r I V ) 


C 

C 

C 

C 


THIS SUBROUTINE LOCATES THE FRAME WITH THE GIVEN TIME 
TIME TAG. IF NO FRAME CAN BE FOUND. A -1 IS RETURNED 




308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 


C AS THE FRAME NUMBER* 

C 

C DECLARE VARIABLES AND DIMENSION ARRAYS* 
C 

INTEGERitc2 IFR<521 ) f IFRD<2) 
EQUIVALENCE < IVALr IFRD< 1 ) ) 

DIMENSION SFRAMEdOOrS) 


OF 


c 

C GET THE TIME REQUESTED IN SECONDS* 

C 

TIME-SFRAME < 1 1 r 1 ) «3600 ♦ 0+SFRAME < 1 1 f 2 ) *60 * 0+SFRAME 
*(II»3)-0*00S 
C 

C SET FRAME NOT READY FLAG* 

C 

NFR»-1 




C 

C GO THROUGH THE RECORDS AND FIND THE TIME. 

C 

K«NREC-1 
DO 100 I«lrK 
IV«I 

READ<IDISK> IFR 
IFRD<1)=IFR<3) 

IFRD<2)»IFR<4) 

C 

C GET THE TIME FOR EACH RECORD* 

C 

ISEC=MAND< IVALf 131071 ) 

IMIL=IFR(5) 

TIMES»FL0AT(ISEC)+*001*FL0AT(IMIL) 

IV»I+1 

READ<IDISK) IFR 
IFRD<1)«IFR<3) 

IFRD<2)=IFR<4) 

ISEC=MAND < IVALf 131071 ) 

IMIL=IFR(5) 

TIMES1«FL0AT<ISEC)+0.001*FL0AT<IMIL) 

NFR=I 

IF(TIME.GE*TIMES .AND* TIME * LT . TIMESl ) GO TO 200 
IF<TIME.GE*TIMES1 .AND* TIME ♦ LT . TIMES ) GO TO 200 
100 CONTINUE 

IF<NFR *EQ* NREC) NFR=-1 
C 

C NO FRAME FOUND* BAIL OUT. 

C 

200 RETURN 
END 

SUBROUTINE RSPRNT<NFRr IDISKf IFR dPRINTr IWORK»NREC»ESNOf 
*IRWD*IPHASEf IV> 

C 

C 

C SUBROUTINE TO PRINT OUT RAW SYMBOLS. 

C 

C DECLARE VARIABLES AND DIMENSION ARRAYS. 

C 

DIMENSION ESN0(8) 

INTEGER*2 IFR<521 ) r IFRD<2) r IW0RK(1025) fMAP 
EQUIVALENCE < IVALr IFRD< 1 ) ) 

DATA IRDFf IDBAR/'SYNC' f 'DB V 


. 368 

3-69 

370 

371 
* 372 

373 

374 
1 375 

376 

377 
, 370 

379 

380 

381 
' 382 

383 

384 
■, 385 

386 

387 

388 
! 389 

1 390 

391 
I 392 
i 393 

394 

395 

396 

397 

398 
f 399 
, 400 

401 
, 402 

403 
‘ 404 

405 
I 406 
1 407 


408 

409 

410 

411 

412 

413 

414 

415 

416 


i 



417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 


C FORMAT STATEMENTS USED BY ROUT-NE. OP POOR QUAUi"/ 

C 

9 FORMAT <lH0r//6Xr 'PRINTOUT OF RAW CHANNEL 'f 
5K' SYMBOLS FOR SYNCHRONIZED DATA FRAME 'r 
!|ci4f//6Xf 'DATA START TIME-' r/lOXr 'HOURS- 't 
5KI2f/10Xr 'MINUTES- ' r 12 »/10Xr 'SECONDS- 'fF6«3) 

10 F0RMAT<lH0r//6Xf 'PRINTOUT OF RAW CHANNEL '» 

)K' SYMBOLS FOR TAPE FRAME (RECORD) '» 14 r//6Xr 'DATA START TIME-' f 
#/10Xr 'HOURS- ' f I2r/10Xf 'MINUTES- 'k 12 f/lOXr 'SECONDS- 'fF6*3) 

11 FORMAT <lH0r/6Xf 'DATA STOP TIME-' »/10Xr 'HOURS- 'rI2f/10Xr 
3|C' MINUTES- ' f 12 f/lOXr 'SECONDS- 'rF6»3> 

12 FORMAT <1H0»/6X» 'RAW CHANNEL SYMBOLS 'f 
♦'(SOFT-DECISIONS IN OCTAL FORMAT)-') 

13 FORMAT (lH0r//6X» 'RAW CHANNEL SYMBOLS 'r 
♦'(HARD-DECISIONS IN BINARY FORMAT)-') 

14 F0RMAT(/3Xr64Ilr 

♦/3X/64Ilf/3X»64Ilr/3Xf64Ilf/3Xf64Ilr/3Xf64Ilr 
♦/3X»64Ilr/3Xr64Il»/3Xf64Ilr/3Xr64Ilf/3X»64Ilr 
♦/3X»64Ilf/3Xf 6411 »/3Xr64Ilf/3X>64Ilf/3Xf 6511) 

15 FORMAT (IHOrSXf 'SYNC WORD LOCATION-' r 15) 

16 FORMAT (lH0f5Xr'SI8NAL-T0-N0ISE RATIO (DB)-' r 017* 10) 

17 F0RMAT(lH0f/5Xf 'SYNCHRONIZED RAW CHANNEL SYMBOLS 'r 
♦'(SOFT-DECISIONS IN OCTAL FORMA''-') 

18 F0RMAT(lH0r//5X» 'SYNCHRONIZED CHANNF.L SYMBOLS 'r 
♦ '(HARD-DECISIONS IN BINARY FOKM' D-') 


C 

C GET START TIME FOR THE FRAME AND PRINT IT OUT* 

C 

IF(NFR*LT.1*0R.NFR*GT.NREC) GO TO 200 
IV-NFR 

READ(IDISK) IFR 
IFRD(1)-IFR(3) 

IFRD(2)=IFR(4) 

ISEC=MAND( IVALf 131071 ) 

IMIL-IFR(S) 

IHOUR-ISEG/3600 
ISEC-ISEC- ( IH0UR^3600 ) 

IMIN-ISEC/60 

ISEC=ISEC-(IMIN^60) 

SECOND-FLOAT ( ISEC ) + * OOl^FLOAT ( IMIL ) 

C 

C GET THE CHANNEL SYMBOLS UNPACKED BEFORE WRITING OVER THE FRAME* 
C THEN GET STOP TIME* 

C 

C FIRST CHECK IF PRINTED RAW SYMBOLS ARE TO BE SYNCHRONIZED* 

C 

IFdRWD *NE* IRDF) GO TO 20 

WRITE ( IPRINT f 9 ) NFR » IHOUR » IMIN r SECOND 

GO TO 30 

20 WRITE( IPRINT f 10) NFRf IHOUR r IMINr SECOND 
CALL UNPCKS(IFRflWORK) 

30 IF(NFR+1*GT*NREC) GO TO 100 
IV-NFR+1 
READ(IDISK) IFR 
IFRD(1)=IFR(3) 

IFRD(2)=IFR(4) 

ISEC=MAND( IVAL»131071 ) 

IMIL-IFR(5) 

IHOUR-ISEC/3600 
ISEC=ISEC-< IH0UR^3600) 

IMIN-ISEC/60 

ISEC=ISEC-(IMIN^60) 

SEC0ND*^'FL0AT(ISEC) + *001^FL0AT(IMIL) 

WRITE( IPRINT r 1 1 ) IHOUR » IMIN r SECOND 



IF<IRWD ^NE. IRDF) GO TO 100 
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1 


I 


t 


I 


1 

i 

1 

1 

1 

i 

i 

I 
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436 
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440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 

451 

452 

453 

454 

455 

456 

457 

458 

459 

460 

461 

462 

463 

464 

465 

466 

467 

468 

469 

470 

471 

472 

473 

474 

475 

476 

477 

478 

479 

480 

481 

482 

483 

484 

485 

486 

487 

488 

489 

490 

491 

492 

493 


C 
C 

C WRITE OUT CHANNEL SYMBOLS AS THEY APPEAR IN A SYNCHRONIZED 
C DATA FRAME, PIRST UN-MAP THE SYMBOLS. 

C 

IFdPHASE .EQ. IDBAR) GO TO 40 

DO 35 I>=1»1024 

M»IUORK(I> 

I1-12-M '^L L . rj 

IF<M .LE. 4> Il»4-M OF OUAyW 

IW0RK<I)»I1 
35 CONTINUE 
GO TO 46 

40 DO 45 1*1 r 1024 

M»IWORK<I) 

Il»M-5 

IF(M .LE. 4) Il»M+3 
IW0RK(I)»I1 

45 CONTINUE 
C 

C WRITE OUT CHANNEL SYMBOLS FOR (SYNC) FRAME. 

C 

46 WRITE(IPRINT»15) IW0RK(1025) 

WRITE(IPRINTrl6) ESNO(l) 

WRITE(IPRINT»17) 

WRITE(IPRINTrl4) (XWORiX<I)r 1 = 1 r 1024) 

DO 50 1=1 f 1024 

IWORK ( I ) = IHARD < MAP < I WORK < I ) ) ) 

50 CONTINUE 

WRITEdPRINTf 18) 

WRITE(IPRINTrl4) (IWORKd). I = ld024) 

RETURN 

C 

C WRITE OUT CHANNEL SYMBOLS FOR TAPE RECORD. 

C 

100 WRITEdPRINTrl6) ESNOd) 

WRITEdPRINTd2) 

WRITEdPRINTf 14) IWORK 

DO 150 1=1 f 1025 

IWORKd) = IHARD(MAPdWORKd) ) ) 

150 CONTINUE 

WRITEdPRINTf 13) 

WRITEdPRINTf 14) IWORK 
200 RETURN 
END 

SUBROUTINE UNPCKdFRf IWORK) 

C 

C THIS SUBROUTINE UNPACKS A FRAME OF CHANNEL SYMBOLS 
C FROM A DATA RECORD. 

C 

C DIMENSION ARRAYS AND DECLARE VARIABLES. 

C 

INTEGER5IC2 IFR (521 ) f IWORK < 1025 ) dSOf ISl f IS2f MAP 
C 

C GO THROUGH THE DATA RECORD AND UNPACK EACH SYMBOL. 

C 

DO 100 I=lf512 
INDEX=I+9 

CALL CHSYMdFRdNDEX) f ISOf ISlf IS2) 

IND=2*I-1 


0 


494 

495 

496 

497 

498 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 


OF POOR QLVjUTY 


IWORK<IND)=MAP<ISO) 

IWORK< IND+1 )«MAP<IS1 ) 

100 CONTINUE 

CALL CHSYH<IFR<521)fIS0rISlrIS2) 
IW0RK(1025>»MAPaS2) 

RETURN 

ENTRY UNPCKS<IFR»IWORK) 

1^0 200 Ialf512 
lNDEX»I+9 

CALL CHSYM < IFR < INDEX ) r ISO r ISl » IS2 > 
INDa2tI-l 
IWORK(IND)aISO 
IWORK< IND+1 )=IS1 
200 CONTINUE 

CALL CHSYM<IFR<521)fIS0fISlrIS2) 

IU0RK<1025)*IS2 

RETURN 

END 




Appendix B 

Fano Algorithm Subroutine Fortran Code 


9 


OF POOR 


1 SUBROUTINE FANO(ESNODBr INr I0UT» ITCTr ACBr NSCx IDFFLG» ITPFLG) 


4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 


C 

C FORTRAN FANO SEQUENTIAL riECODER FOR OCTAL CHANNEL — 

C 

C THIS PROGRAM SEQUENTIALLY DECODES FRAMES OF 1024 CHANNEL 

C SYMBOLS USING SOFT-DECISIONS ON THE DEMODULATED CHANNEL 

C OUTPUT. EACH FRAME CORRESPONDS TO 512 INFORMATION BITS OF 

C WHICH THE LAST 24 IS A KNOWN FRAME SYNC WORD WITH THE FOLLOWING 
C pattern: 111110001100010101001001. this known tail sequence 

C IS USED BY THE FANO ALGORITHM IN THE DECODING PROCESS. 

C 

C IN = INPUT ARRAY OF 1024 SOFT-DECISIONS THAT ARE INPUT 
C TO THE FANO DECODER 

C lOUT = OUTPUT ARRAY OF 512 DECODED INFORMATION BITS 
C ITCT = ON ENTRANCE TO ROUTINEf IS THE MAXIMUM NUMBER 
C OF COMPUTATIONS ALLOWED PER FRAME 

C ACB = AVERAGE NUMBER OF COMPUTATIONS PER INFORMATION BIT 

C DECLARED AS BEING DECODED. 

C NSC = NUMBER OF CORRECTED CHANNEL SYMBOLS 

C < =0 WHEN FRAME IS DELETED) 

C IDFFLG = DELETED FRAME FLAG 

C 0 = FRAME SUCCESSFULLY DECODED 

C 1 = FRAME DELETED - MAXIMUM NUMBER OF COMPUTATIONS 

C ALLOWED PER FRAME HAS BEEN EXCEEDED 

C 2 = INVALID DECODER PARAMETER 

DIMENSION LK600> »IREG<600) fLMDA<600) r IPO<600) rIPK600) 

* fQL(7) fTP<8) fD<8) .'QV<64) » MOO (64) fM 01<64) i-ITAIL<36) 

* fINHD(1500) 

INTEGER*2 IN(1024) f I0UT(512) 

C 

C ESNODB = SIGNAL ENERGY TO NOISE DENSITY RATIO IN DB 

C QSP = QUANTIZATION SPACING 

C (NORMALIZED TO STANDARD DEVIATIONS OF THE NOISE) 

C IFL = TOTAL FRAME LENGTH 

C ITL = LENGTH OF FIXED TAIL SEQUENCE 

C ITRIMA = MAXIMUM TRIALS PER BLOCK 

C (SET EQUAL TO ITCT INITIALLY) 

C ITAPOArlTAPlA - REGISTER TAP MASKSr LEFT JUSTIFIED 

C BIAS = METRIC BIAS PER BRANCH 

C TO = THRESHOLD SPACING 

C SCALE = SCALE FACTOR FOR INTEGER COMPUTATIONS 

ITRIMA=ITCT 
DATA ITP /'TAPE'/ 

DATA QSP /0.5/ 

DATA ITAPOAf ITAPIA /ZADD6F7DDr ZEDD6F7DD/ 

DATA BIAS /l.O/ 

DATA TO /3.0/ 

DATA SCALE /lOOO.O/ 

DATA IFL /512/ 

DATA ITL /24/ 

DATA ITAIL(l) f ITAIL(2) rITAIL(3) i-ITAIL(4) /Irltlfl/ 

DATA ITAIL(5) f ITAIL(6) »ITAIL(7) f ITAIL(8) /l.OfOrO/ 

DATA ITAIL(9) f ITAIL(IO) f ITAIL(ll) 7lTAIL(12) /Irl^OrO/ 

DATA ITAIL(13) rITAIL(14) ,ITAIL(15) rITAIL(16) /OflrOfl/ 

DATA ITAIL(17)fITAIL(18).ITAIL(19)fITAIL(20) /O.lyO.O/ 

DATA ITAIL(21) f ITAIL(22) rITAIL(23) .ITAIL(24) /1.0..0fl/ 

DATA RAL0G2 /I. 442695/ 


58 C 

59 ISIGN=-1073741824 

60 ISIGN=ISIGN+ISIGN 

61 MAX=2*IFL 


• 

62 

c 


, • *1 

• <• 



63 

c 


PROGRAM INITIALIZATION 

, 


64 

c 


INITIALIZE CHOICE ARRAY 



65 

c 


LI = 0 FOR BEST CHOICE <0 ON EOEN NODE TIEf 1 ON OBD) 



66 

c 


1 FOR WORST CHOICE (1 ON EVEN NODE TIE> 0 ON ODD) 



67 

c 


2 FOR KNOWN BIT IN TAIL SEQUENCE 


1 

68 



DO 32 I=lfIFL 



69 


32 

LI(I)=0 



70 



IFLM1=IFL-1 



71 



IFLP1=IFL+1 


1 

72 



ITB=IFLP1-ITL 

* 


73 



DO 33 I=ITBrIFLPl OF POOR QUAUTY 



74 


33 

LKI)=2 


1 

75 



ITO=TO*SCALE 



76 



ESN0=>EXP < ESNODB/10 ♦ OJKALOG < 10 ♦ 0 ) ) 



77 



TH=SQRT<2.0*ESN0) 


. 

78 

c 





79 

c 


CALCULATE CHANNEL TRANSITION PROBABILITIES 



80 



DO 34 I=lf7 



81 


34 

QL < I ) =TH-FLOAT ( 4-1 ) litaSP 


1 

82 



QV(1)=0*5*<2.0-ERFC<0«7071068J»CQL(1) ) ) 


j 

83 



TP(1)=QV(1) 


1 

84 



DO 35 1=2 r 7 


1 

85 



Q0(I)=0«5*<2.0-ERFC(0.7071068NcQL(I) ) ) 


86 


35 

TP<I)=QV<I)-Q0(I-1> 


I 

87 



TP(8)=1*0-QV(7) 



88 



PE=1,0-Q0<4> 


• 

89 

c 





90 

c 


CALCULATE BRANCH METRICS 



91 



DO 36 l=lf3 


t 

92 



IM9-9-I 



93 


36 

D(I)=AL0Q(TP(I)*2»0/(TP(I)+TP(IM9) ) ))tcRAL0G2 


t 

94 



DO 37 1=1 »8 



95 



DO 37 J=lf8 



96 



MET= < D ( I ) +D < J > -BI AS ) >lcSCALE 



97 



K=(I-1)*8+J 



98 



M00<K)=MET 



99 



K=(I-l)!|c8+9-J 


i 

100 


37 

M01(K)=MET 


1 

101 

c 





102 

c 


OBTAIN HARD-DECISIONS FROM RECEIVED CHANNEL SYMBOLS 



103 



JO AX J=lfMAX 



104 



IF(IN<J) .GE. 5) GO TO 40 



105 



INHD< J)=0 



x06 



GO TO 41 



107 


40 

INHD< J)=l 


- 

108 


41 

CONTINUE 



109 

c 





110 

c 


INDEX RECEIVED SYMBOLS (SOFT-DECISIONS) INTO BRANCH METRIC ARRAY 


1 

ill 

c 


EACH PAIR OF SYMBOLS DEFINES AN ARRAY ELEMENT 


- 

112 



K=0 



113 



DO 50 I=l»MAXf2 


t 

114 



INDEX=8# ( IN ( I ) -1 ) +IN ( I+l ) 



115 



K=K+1 



116 



IN(K)=INDEX 



117 


SO 

CONTINUE 


1 

ilS 

c 





11? 

c 


INITIALIZE DECODER BEFORE START OF EACH FRAME 



120 



IP0A=0 


1 

121 



IP1A=0 



122 



ITCT=0 



123 



IT=0 


• 

124 



L=0 

i-l 



' n - 



- 4,^ 



125 



IFLAG=0 ORIGINAL PAGE S3 


126 



N=1 OF POOR QUALITY 

•1 

127 



LK1)=0 


12S 

C 




129 

C 


DECODE IFL NODES 


130 


80 

IF<N-IFL) 81f81f310 

1 

131 


81 

IND=IN<N) 


132 



INDI=65-IND 


133 



IF<IP0A) 135fl31fl31 

J 

134 


131 

IF(IPIA) 134»133»133 


135 


133 

M0=M00<IND) 


136 



M1=M00(INDI) 


137 



GO TO 139 

1 

138 


134 

MO=MOKIND) 


139 



M1=M0KINDI) 


140 



GO TO 139 

1 

141 


135 

IFCPIA) 137fl36»136 


142 


136 

M0=M01(INDI) 


143 



M1=M0KIND) 

1 

; 

144 



GO TO 139 

145 


137 

M0=M00<INDI) 

i 

146 



M1=M00<IND) 


147 

C 



1 

148 

C 


CHECK FOR KNOWN BIT 

i 

149 


139 

IF<LKN)-2)140f 165»999 


150 


140 

IF<M0-M1) 145fl41»149 

1 

151 


141 

IF(N-2*<N/2) ) 999»149rl45 

J 

152 


145 

IF(LKN)) 999..160fl70 

j 

153 


149 

IF<LI(N)) 999fl70fl60 


154 

C 



1 

155 

C 


TRIAL DECISION IS A ONE 

a 

156 


160 

LMDA(N)=M1 


157 



IREG(N)=1 

1 

158 



LT=L+M1 


159 



GO TO 210 


160 

C 



1 

161 

C 


KNOWN TAIL SEQUENCE f DO NOT TIGHTEN THRESHOLD 

1 

162 


165 

IARG=N-ITB+i 

« 

163 



IF(ITAILdARG) ) 999 f200f201 


164 


200 

LMDA<N)=MO 

1 

165 



IREG<N)=0 

1 

166 



LT=L+MO 

167 



IF(LT-IT) 270f166f166 


168 


201 

LMDA(N)=M1 


169 



IREG(N)=1 


170 



LT=L+M1 


171 



IF(LT-IT) 270fl66Fl66 


172 


166 

IF(IFLAG) 999 r 230 f 250 

a 

173 

C 




174 

C 


TRIAL DECISION IS A ZERO 


175 


170 

LMDA<N)=MO 


176 



IREG(N)=0 

1 

177 



LT=L+MO 


. 178 


210 

IF(LT-IT) 270f219f219 


179 


219 

IFdFLAG; 999 f 220 » 250 

* 

ISO 

c 




131 

c 


CHECK IF THRESHOLD CAN BE TIGHTENED 


182 

c 


UPDATE TRIALS COUNT 

V 

183 


220 

IF(LT-<IT+ITO) ) 230f221f221 


184 


221 

IT=IT+ITO 


185 


230 

L=LT 


186 



irCT=ITCT+l 


187 


IF(IREG<N)) 999r232»231 

188 

231 

IPlAaMXOR(IPlAiITAPlA) 

189 

232 

IF(IPIA) 234 » 233 r 233 

190 

233 

IPlA=MOVEL<IPlAfl> 

191 


IP1(N)=0 

192 


GO TO 235 

193 

234 

IP1A=M0VEL<IP1A»1) 

194 


IP1(N)=ISIGN 

195 

235 

IF(IREG<N)) 999»237»23o 

196 

236 

IPOA=MXOR < IPOA f ITAPOA ) 

197 

237 

IFdPOA) 239»238f238 

198 

238 

IP0A=M0VEL<IP0A>1) 

199 


IP0<N)=0 

200 


GO TO 242 

201 

239 

IP0A=M0VEL<IP0A»1) 

202 


IPO(N)=ISIGN 

203 

C 


204 

C 

PROCEED TO NEXT NODE 


OF POOR QUAin’^ 


205 

206 

207 

208 

209 

210 
211 
91 9 

«>• X Am 

213 


242 N=N+1 

IF<LKN)-2) 243f80»999 

243 LI(N)=0 
GO TO 80 

C 

C CHECK IF RUNNING THRESHOLD IS VIOLATED 

250 IF<LT-<IT+ITO) ) 260»230»230 
260 IFLAG=0 
GO TO 230 


214 

215 

216 
217 


C 

C SEARCH OTHER LIKELY PATHS - STEP DECODER BACKWARDS 

270 IFLAG=1 
280 ITCT=ITCT+1 


218 IFdTCT-ITRIMA) 282»440»440 

219 282 IF<N-1) 999f300r284 


220 

221 

999 


284 LB=L-LMDA(N-1) 
IF<LB-IT) 300f285f285 

285 N=N-1 


223 

224 
99S 

226 

/^27 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 
233 

239 

240 

241 

242 

243 

244 

245 

246 

247 


L=LB 

IF<IREG<N)) 999f286f287 
286 IPlA=MXOR<MOVER<IPlAf 1) f IPl(N) ) 
IP0A=MX0R<M0VER(IP0Af 1) f IP0(N> > 

GO TO 296 

237 IPlA=MXOR(MXOR<MOVER<IPlArl) f IPKN) ) rITAPlA) 
IPOA=MXOR <MXOR< MOVER < IP0A» 1) » IPO (N) ) rITAPOA) 

296 IF<LKN)) 999»297f280 

297 LKN) = 1 
GO TO 81 

C 

C RELAX RUNNING THRESHOLD VALUE 

300 IT=IT-ITO 
IF<LKN)-2) 301»80f999 

301 LKN)=0 
GO TO SO 


C 

C 

310 


C 

C 

C 

C 

c 

c 


IFL NODES HAVE BEEN DECODED OR TRIALS LIMIT EXCEEDED 
N=N-1 

ITCT = ON EXIT FROM ROUTINE» NUMBER OF FORWARD AND 
REVERSE MOVES MADE DURING DECODING PROCESS 
IT = FINAL THRESHOLD VALUE (SCALED) 

L = FINAL PATH METRIC VALUE (SCALED) 

LMDA = PATH METRIC (SCALED) 

IREG(N) = DECODED DATA 


G4 w w w r j ro r o r o fo ro r j ro r-j ro r o r j r j tj r j r j ro ro r o r o ro r j ro ro ro ro ro r o ro r j ro to to ro to to to to to to to to to to to to to to to to to to 
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ORIGSS'^AL 
OF POOs^ QUA'i 


C CONVERT BRANCH SYMBOL ARRAYS TO 0»1 FORMAT 
DO -416 Jsl,N 
IF<IP0(J)) A10rA12f999 
410 IP0<J)=1 

412 IF<IPKJ)) 414M16r999 
414 IPKJ) = 1 
416 CONTINUE 

c 

C CALCULATE THE AVERAGE NUMBER OF COMPUTATIONS PER 
C DECODED INFORMATION BIT 

ACB=FLOAT<ITCT)/FLOAT<IFL) 

C 

C DETERMINE THE NUMBER OF CORRECTED CHANNEL SYMBOLS 
NSC«0 
K=1 

DO 424 J=lfN 

IF(INHD<K)-IPO< J) ) 420r421f420 

420 NSC=NSC+1 

421 IF<INHD(K+1)-IP1( J) ) 422^423^422 

422 NSC=NSC+1 

423 K=K+2 

424 CONTINUE 

DO 430 J=1fN 
IOUT(J)=IREG< J) 

430 CONTINUE 
C 

C NORMAL EXIT FROM DECODER 

IDFFLG=0 
RETURN 

440 CONTINUE 
C 

C DELETE FRAME - MAXIMUM NUMBER OF COMPUTATIONS ALLOWED 

C HAS BEEN EXCEEDED 

N=N-1 
IDFFLG=1 
NSC=0 
C 

C CALCULATE AVERAGE NUMBER OF COMPUTATIONS PER 

C BIT DECLARED AS BEING DECODED 

ACB=FLOAT(ITCT)/FLOAT(N) 

IFdTPFLG .EQ* ITP) GO TO 450 
DO 441 J=lfN 
I0UT( J)=IRE0(J) 

441 CONTINUE 
M=N+1 

no 442 J«MrS12 
I0UT< J)=9 

442 CONTINUE 
RETURN 

450 ILT=IFL-ITL 

DO 451 J=lfILT 
I0UT( J)=0 

451 CONTINUE 

DO 452 J=lfITL 
I0UT< J+488)=ITAIL( J) 

452 -CONTINUE 
RETURN 


306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 


C 

C DELETE FRAME - INVALID DECODER PARAMETER WAS 
C ENCOUNTERED DURING EXECUTION 
999 IDFFLG=2 
NSC=0 
ACB=0»0 
ILT=IFL-ITL 
DO 500 J=lrILT 
I0UT<J)=0 
500 CONTINUE 

DO 510 J=lrITL 
I0UT(J+488)=ITAIL<J) 

510 CONTINUE 
RETURN 
END 


ORJGJrjAL as 

OF POOR OU/V„(TY 


Appendix C 

Utility Subroutine Assembler Code 




1 





2 





3 

BEGN 

CSECT 

READONLY 


« A 


ENTRY 

MXOR f NOVEL f MOVER r MAND f MLSL f MLSR 

5 


ENTRY 

M0R»M0R2^MAND2 

fMLSL2fMLSR2fIFLIP 

6 

MXOR 

B 

5+5<15> 

ROUTINE TO EXCLUSIVE OR THE 

, 7 


DC 

X'5' 

CONTENTS OF TWO FULLWORDS# 

a 

■ 

DC 

CLS'MXOR" 


9 


SAVE 

<14f 12) 


10 


USING 

MXORrlS 


^ 11 


L 

2f0<l) 


12 


L 

0f0<2) 

OF POOR QUALH’Y 

13 
. 14 


1 . 

(- 

1>4<1) 

IfO(l) 

15 


XR 

Ofl 


16 


LA 

ISrO 


• 17 


LM 

2rl2f28(13) 


18 


MVI 

12(13) rX'FF' 


19 


BR 

14 


20 

* 




1 21 





i 22 

MOOEL 

B 

5+5(15) 

ROUTINE TO LOGICAL SHIFT LEFT ONE 

23 


DC 

X'S" 

POSITION THE CONTENTS OF A FULLWORD 

, 24 


DC 

CL5' HOVEL' 


1 25 


SAVE 

.14»12) 


J 26 


USING 

NOVEL r 15 


27 


L 

IrO(l) 


f 28 


L 

0»0(1) 


1 29 


SLL 

Orl 


30 


LA 

15f0 


. 31 


LM 

2f 12^28(13) 


32 


MVI 

12(13) fX'FF' 


' 33 


BR 

14 


34 





’ 35 

* 




1 36 

MOOER 

B 

5+5(15) 

ROUTINE TO LOGICAL SHIFT RIGHT ONE 

37 


DC 

X'5' 

POSITION THE CONTENTS OF A FULLWORD 

.] 38 


DC 

CL5' MOVER' 


j 39 


SAVE 

<14f 12) 


^ 40 


USING 

MOVERrlS 


J 

1 42 


L 

IrO(l) 



L 

0»0(1) 


1 43 


SRL 

Ofl 

• 

44 


LA 

15f0 


45 


LM 

2rl2r28(13) 


46 


MVI 

12(13) fX'FF' 


■' 47 


BR 

14 


48 

* 




49 

t 




• 50 

HAND 

B 

5+5(15) 

ROUTINE TO AND THE CONTENTS 

* 51 


DC 

X'5' 

OF TWO FULLWORDS* 

52 


DC 

CLS'MAND' 


4 53 


SAVE 

(14f 12) 


54 


USING 

MANDrlS 


55 


L 

2f0(l) 


56 


L 

0r0(2) 


■ 57 


L 

lf4(l) 


58 


L 

IfO(l) 


59 


NR 

Ofl 


• 60 


LA 

15f0 


61 


LM 

2f 12f2S(13) 


62 


MVI 

12(13) fX'FF' 


63 

64 

* 

BR 

14 



Vf’ ■ •• • 4 ■ ■ . 



■y „.' ^ -■■■■ " ■■ xY ' 


65 

* 

. 


I*- 

' 66 

HLSL 

B 

5+5<15) 

ROUTINE TO LOGICAL SHIFT LEFT THE 

67 


DC 

X'5' 

CONTENTS OF A FULLWORD A SPECIFIED 

66 


DC 

CLS'MLSL' 

NUMBER OF POSITIONS. | 

* 69 


SAVE 

<14rl2> 

ii 

70 


USING 

MLSLf 15 

ii 

ii 

71 


BALR 

9r0 

Ii 

, ^2 


L 

2f0(l) 


73 


L 

0f0<2) 


74 


L 

lf4<l) 


75 
’ 76 


L 

LA 

IrO(l) 

2r0 

ORIGIMAI 53 ■ 

77 

LOOPL 

SLL 

Orl 

OF POOR QUALITY ' i 

. 78 


LA 

2rl<2> 


. 79 


CR 

lr2 

, 

80 


BNE 

LOOPL 


81 


LA 

15f0 


. 82 


LM 

2rl2»2Sa3) 


' 83 


MVI 

12(13)fX'FP' 

• 

84 


BR 

14 


85 
1 86 
1 97 
88 

* 

MLSR 

B 

5+5(15) 

ROUTINE TO LOGICAL SHIFT RIGHT THE 


DC 

X'5' 

CONTENTS OF A FULLWORD A SPECIFIED 

89 
' 90 


DC 

CLS'MLSR' 

NUMBER OF POSITIONS. 


SAVE 

<14f 12) 


91 


USING 

MLSRf 15 


92 


L 

2»0(1) 


• 93 


L 

0?C<2) 


94 


L 

lf4(l) 


95 


L 

IfO(l) 


. 96 


LA 

2f0 


97 

LOOPR 

SRL 

Qfl 


1 98 ■ 


LA 

2rl<2) 


99 


CR 

lr2 


100 


BNE 

LOOPR 


. 101 


LA 

15f0 


102 


LM 

2f 12 p2S(13) 


103 


MVI 

12(13) fX'FF' 


104 


BR 

14 


105 

* 




106 

* 




107 

MOR 

B 

3+5(15) 

ROUTINE TO INCLUSIVE OR THE 

108 


DC 

X'3' 

CONTENTS OF TWO FULLWORDS. 

109 


DC 

CL3'M0R' 


110 


SAVE 

(14»12) 


111 


USING 

MORf 15 


•’ 112 


L 

2»0(1) 


113 


L 

OfO<2) 


114 


L 

lf4(l) 


' 115 


L 

IrO(l) 


. 116 


OR 

Of 1 

INCLUSIVE OR 

117 


LA 

15f0 


. 118 


LM 

2f 12f28(13) 


119 


MVI 

12(13) fX'FF' 


120 


BR 

14 


121 





' 122 

* 




123 

M0R2 

B 

5+5(15) 

ROUTINE TO INCLUSIVE OR THE CONTENTS 

124 


DC 

X'5' 

OF A HALFWORD AND A FULLWORD. 

- 125 


DC 

CL5'M0R2' 


126 


SAVE 

(14fl2) 


127 


USING 

M0R2f 15 


128 


L 

2f0(l) 

GET ADDRESS OF FIRST OPERAND 

129 


LH 

0f0(2) 

LOAD REGO WITH DATA IN HALFWORD | 



30 


N 

0f*X'0000FFFF 

131 


L 

1»4<1) 

,132 


L 

lfO<l> 

.33 


OR 

Orl 

x34 


LA 

15r0 

135 


LM 

2rl2f28(13) 

* 36 


MVI 

12(13) rX'FF' 

37 


BR 

14 

138 




i'39 

* 



.40 

MAND2 

B 

5+5(15) 

141 


DC 

X'5' 

,142 


DC 

CL5'MAND2" 

.43 


SAVE 

(14»12) 

-44 


USING 

MAN02rlS 

. 145 


L 

2»0(1) 

' .46 


LH 

0r0(2) 

.47 


L 

lf4(l> 

148 


LH 

1.0(1) 

,1.49 


NR 

Of 1 

, .SO 


N 

Of*X"OOOOFFFF 

*151 


LA 

15f0 

152 


LM 

2f 12.28(13) 

! L53 


MVI 

12(13) .X^FF' 

1 L54 


BR 

14 

155 

* 



1^56 

♦ 



.57 

MLSL2 

B 

5+5(15) 

'l58 


DC 

X'S' 

.159 


DC 

CL5'MLSL2' 

' .60 


SAVE 

(14.12) 

<_61 ■ 


USING 

MLSL2.15 

162 


L 

2.0(1) 

' .63 


LH 

0.0(2) 

.64 


N 

0.=X'0000FFFF 

165 


L 

1.4(1) 

166 


L 

1.0(1) 

.67 


LA 

2.0 

‘i68 

L00PL2 

SLL 

0.1 

169 


LA 

2.1(2) 

1 170 


CR 

1.2 

( 171 


BNE 

LOOPL2 

172 


LA 

15.0 

173 


LM 

2.12.28(13) 

174 


MVI 

12(13) .X'FF' 

.• 175 


BR 

14 

176 

* 



177 

* 



■ 178 

MLSR2 

B 

5+5(15) 

. 179 


DC 

X'5' 

180 


DC 

CL5'MLSR2' 

, 131 


SAVE 

(14.12) 

182 


USING 

MLSR2.15 

183 


L 

2.0(1) 

184 


LH 

0.0(2) 

“ 185 


N 

0.=X'OOOOFFFF 

186 


L 

1.4(1) 

187 


L 

1.0(1) 

. 188 


LA 

2.0 

189 

L00PR2 

SRL 

0.1 


GET ADDRESS OF SECOND OPERAND 
LOAD REGl WITH DATA IN FULLWORD 
INCLUSIVE OR 


ROUTINE TO AND THE CONTENTS 
OF TWO HALFWORDS* 


. ’Ji ■ 

OF POOr^ 


ROUTINE TO LOGICAL SHIFT LEFT THE 
CONTENTS OF A HALFWORD A SPECIFIED 
NUMBER OF POSITIONS. 


GET ADDRESS OF FIRST OPERAND 
LOAD REGO WITH DATA IN HALFWORD 

GET ADDRESS OF SECOND OPERAND 
LOAD REGl WITH SHIFT DATA 

SHIFT LEFT ONE POSTION 
INCREMENT LOOP COUNTER 
CHECK IF SHIFTED ENOUGH 


ROUTINE TO LOGICAL SHIFT RIGHT THE 
CONTENTS OF A HALFWORD A SPECIFIED 
NUMBER OF POSITIONS. 


■i . : , 4 


0R!«?3tWAL‘ 

OF POOR QUAUHC 




190 


LA 

2fl<2> 

191 


CR 

lf2 

♦ 192 


BNE 

L00PR2 

193 


LA 

ISrO 

194 


LM 

2rl2f28(13> 

, 195 


MVI 

12(13)fX'FF 

196 


BR 

14 

197 

* 



198 

* 



• 199 

IFLIP 

B 

5+5<15) 

200 


DC 

X'S" 

201 


DC 

CL5' IFLIP" 

202 


SAUE 

<14rl2) 

203 


USING 

IFLIPrlS 

204 


L 

2f0<l) 

. 205 


LH 

1»0<2> 

' 206 


L 

2fBMl 

\ 207 


NR 

2rl 

208 


SLL 

2f 11 

f 209 


L 

3fBM2 

3 210 


NR 

3f 1 

211 


SRL 

3rll 



L 

4rCBMl 

213 


NR 

if4 

214 


XR 

lf3 

215 


L 

4fCBM2 

216 


NR 

lf4 

. 217 


XR 

lr2 

218 


LR 

Ofl 

, 219 


LA 

15f0 

220 


LM 

2f 12^28(13) 

OOI 


MVI 

12(13) »X"FF 

222 


BR 

14 

} 223 

BMl 

DC 

X" 00000007" 

1 224 

BM2 

DC 

X" 00003800" 

225 

CBMl 

DC 

X"FFFFFFF8" 

226 

CBM2 

DC 

X"FFFFC7FF" 

227 


END 



ROUTINE TO SWITCH -THE POSITIONS OF 
THE 3-BIT SOFT-DECISIONS PACKED 
INTO A HALFWORD. 


GET ADDRESS OF OPERAND 
GET HALFWORD OPERAND 
GET MASK FOR BITS 14-16 
REG2«BITS 14-16 
SHIFT BITS TO NEW POSITION 
GET MASK FOR BITS 3-5 
REG3»BITS 3-5 

SHIFT BITS TO NEW POSITION 

AND OUT OLD BITS 14-16 
OR IN'nEW bits 14-16 

AND OUT OLD BITS 3-5 
OR IN NEW BITS 3-5 


I 




OF POOti 


tn 

n 



228 

* 




t 

229 

* SUBROUTINE 

TO PACK DECODED 

BITS - CALL DPACK(INARR.IOUTAR) 


230 

BRACK 

CSECT READONLY 



231 


B 

5+5<15) 


, 

232 


DC 

X'5" 



233 


DC 

CLS'DPACK' 



234 


SAME 

(14f 12) 



235 


USING DPACKrlS 


f 

236 

* GET 

INPUT-! 

OUTPUT ARRAY ADDRESSES 


237 


L 

OfO<l) 

GET ADDRESS OF INARR 


238 


L 

lf4<l) 

GET ADDRESS OF lOUTAR 

1 

239 


LA 

3r0 

SET REG3-0 

1 

240 


LR- 

4f3 

SET REG4-0 

1 

241 

* INITIALIZE 

OUTPUT ARRAY TO 

ZERO 

1 

242 


LR 

5f 1 

SET REG5«ADDR OF lOUTAR 

243 

LOOPZ 

STH 

3f0<5) 

STORE 0 IN HALFWORD LOCATION 

1 

244 


LA 

5f2<5) 

GET ADDRESS OF NEXT LOCATION 


245 


LA 

4f 1<4> 

INCREMENT LOOP COUNTER 

1 

246 


C 

4f»F'32' 

TEST IF FINISHED 

i 

247 


BNE 

LOOPZ 



248 

JK GET 

DATA FROM INPUT ARRAY 

AND PACK INTO HALFWORDS 

1 

i 

249 


LR 

5.0 

SET REGSsADDR OF INARR 

250 


LR 

7.3 

INITIALIZE WORD LOOP COUNTER 

251 

LOOPU 

LR 

4.3 

INITIALIZE BIT LOOP COUNTER 


252 


LA 

9.1 

SET REG9«1 


253 


LR 

8r3 

INITIALIZE PACKING REGISTER 

i 

254 

LDOPB 

SLL 

8.1 

SHIFT PACKING REGISTER LEFT 


255 


LH. 

6.0(5) 

LOAD INPUT DATA 

I 

256 


CR 

3.6 

CHECK IF B.TT IS A ZERO 


257 


BE 

CONT 


1 

258 


XR 

8.9 

OR 1 TO REGS 

1 

259 

CONT 

LA 

5.2(5) 

GET ADDRESS OF NEXT HALFWORD 

260 


LA 

4.1(4) 

INCREMENT BIT LOOP COUNTER 

1 

261 


C 

4.=F'16' 

TEST IF C0UNTER=16 


262 


BE 

OUTW 

THIS WORD IS PACKED 

1 

263 


B 

LOOPS 

CONTINUE PACKING THIS WORD 

1 

264 

OUTW 

STH 

8.0(1) 

STORE PACKED DATA IN OUTPUT 

265 


LA 

1.2(1) 

GET NEXT OUTPUT WORD ADDRESS 

1 

266 


LA 

7.1(7) 

INCREMENT WORD LOOP COUNTER 


267 


C 

7.»F'32' 

TEST IF C0UNTER=32 

* 

268 


BNE 

LOOPW 

START PACKING NEXT WORD 


269 

* RETURN TO 

CALLING PROGRAM 



270 


LA 

15.0 


1 

271 


LM 

2.12.28(13) 



272 


MMI 

12(13).X'FF' 



273 


BR 

14 



274 


END 




OK!^£MAL PA'SE B 
OF POOR QUALtTY 


275 ^ 

2-^ * SUBROUTINE TO UNPACK SOFT-DECISIONS - CALL CHSYM< IVALr ISO r ISl » IS2) 



CHSYM 

CSECT 

READONLY 





B 

5+5<15) 





DC 

X'5' 



23C 


DC 

CL5'CHa,.M/ 



281 


SAVE 

(14^12) 



282 


USING 

CHSYM f 15 



283 

% GET 

INPUT-OUTPUT PARAMETER ADDRESSES 


284 


L 

2r0<l) 

GET ADDRESS OF 

IVAL 

285 


LH 

0f0<2) 

GET CONTENTS OF 

IVAL 

286 


L • 

2r4<l> 

GET ADDRESS OF 

ISO 

287 


L 

3r8<l) 

GET ADDRESS OF 

ISl 

288 


L 

4f 12(1) 

GET ADDRESS OF 

IS2 

289 

% MASK 

UNWANTED BITS AND SAVE 3-BIT CHANNEL 

. SYMBOLS 

290 


L 

5rMK0 



291 


NR 

5f0 

SAVE ONLY BITS 

3-5 

292 


SRL 

5f 11 

LEFT JUSTIFY 


293 


STH 

5r0<2) 

SAVE IN ISO 

• 

294 


L 

5»MK1 



295 


NR 

5f0 

SAVE ONLY BITS 

6-8 

296 


SRL 

5f8 

LEFT JUSTIFY 


297 


STH 

5f0<3) 

SAVE IN ISl 


298 


L 

5fMK2 



299 


NR 

5f0 

SAVE ONLY BITS 

14-16 

300 


STH 

5f0(4) 

SAVE IN IS2 


301 

% RETURN TO CALLING PROGRAM 



302 


LA 

15f0 



303 


LM “ 

'2rl2r28<13) 



304 


MVI 

12(13) fX'FF' 



305 


BR 

14 



306 


DC 

BL2'0000' 



307 

MKO 

DC 

X' 00003800' 



308 

MKl 

DC 

X' 00000700' 



309 

MK2 

DC 

X' 00000007' 



310 


END 






c 


Appendix D 
Technical Memos from 


NASA Ames and JFL 


JET PROPULSION LABORATORY 


OP POOR 




R. B. MILLER 

IirrEROFFICE MEMORATO « « 
#3384-77-037 HUia ^ y |377 

25 August 1977 


TO: R. B. Miller . 

FROM; J. H. Milcher^/^* 

SUBJECT: Recording Formats for the Pioneer-Venus Reverse Playback Program 


The Telemetry Processor Assembly has available, as inputs from the Symbol Synchronizer 
Assembly (SSA), three formats for symbol data. Each of the three formats have their 
distinctive use in the Telemetry Processor Assembly. However,, in order' to best 
serve the needs of the Pioneer Venus Proj'ect in its decoding of the Reverse Playback 
data, it is proposed that the project select the format, from the three available, that 
best fits their needs. . ■ 


The three formats available are as follows: 
1. UNSYNCHRONIZED 



R. B. Miller 


-2 


•I0.MJ3384-77-037 
24 August 1977 



where Sq, S^, ... are the quantized symbols with weighting as shown in 

Table 1. , * • 



where Sq ... Sn are the sign bits and Mq ... Hn are the 7-bit magnitudes 
associated with the sign bits. 
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f 

, The only remaining unanswered questions pertain to: (1) the desired record 

length in bits (bytes); (2) what additional data, if any, is required. 





Table 1 
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OET PROPULSION LABORATORY 


UPt 

INTEROrnCE MEflORANOUM 
#3384-78-039 

8 August, 1978 


TO: 0, H. 

FROM: «R. L. Murray 

SUBJECT; ' Digital Tapa Format for the Pioneer Venus Reverse PTayback Process. 


The tapes produced’ by the Pioneer Venus Reverse Playback Process v/ill be 
standard 9 track digital tapes, recorded at 1600 bpi, using phase encoding. 

Each tape record will consist of 521 16-bit v/ords (see attached figure). An 
End of File (EOF) mark will be written after the last record on each data tape. 
A detailed description of the tape format follows. 

1. Word 1 contains the tape record number (0,1,2, etc. through 65,536). 

2. Vlord 2 contains the following playback equipment lock status 
information: 

^ Bit 10 » Receiver Lock Status. 

0 => In Lock 

1 = Out of Lock 

Bit 11 = SDA Lock Status 

0 » In Lock 

1 a Out of Lock 

Bit 12 a SSA Lock Status 

0 = In Lock 

1 a Out of Lock 

3. Bits 6 through 14 of word 3 contain the playback day of year 
in binary. 

4. Bit 15 of word 3 and bits 0 through 15 of word 4 cohtain the 
time of day of olayback in binary seconds (B15 of word 3 is the 
MSB and B15 of word 4 is the LSD). 

5. Bits 6 through IS of viord 5 contain the milliseconds of second of 
playback in binary. 
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6. ’ Word 6 contains the follov/ing playback equipment configuration 

information: 

Bits 0 - 3 ■ The receiver number (1-7) 

Bits 4-7 =• Tlie SDA number (1-8) 

* Bits 8 -11 = The SSA number (1 or 2) 

Bits 12-15 « The TPA number (1 or 2) 

7. Word 7 contains the receiver AGC reading in tvra's complement 
binary representation with the binary point located between 
bits 8 and 9. 

8. Word a contains the SSA Signal-to-»Moise Ratio in dbm in 
two's complement binary representation with the binary point 
located betv;een bits 8 and 9. 

9. Vtord 9 contains zeroes. 

10. VJords 10 thi=Qugh 521 contain quantized symbol values associated 
with 512 telemetry data bits (1024 symbols). So is the first 
symbol received by the digital recording program, S- the second, 
etc. 
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E. C. Gatz 
R. B. Miller 
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8 August, 1973 


PIONEER VENUS RF/ERSE PLAYBACK 
DIGITAL TAPE FORMAT 


Vtords / Bits 


Q- 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

RECORD NUMBER 


LOCK STATUS 


3 0 0 0 0 0 0 


DAY OF YEAR 


SECONDS OF DAY 


S 0 0 0 0 0 0 


MILLISECONDS 


RECEIVER AGC 
SSA SNR • 


SPARE 


S2 S3 


( 512 WORDS) 


GO S1023 S0124 


O is • «> . » 
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Enclosure 1 

P-V REVERSE PLAYBACK DECODING 
OUTPUT TAPE FORMAT 


The Output tapes shall be In the standanJ Pioneer Intermediate Data 
Record (IDR) format. These tapes will be used in the Pioneer .Venus Data 
Records Processing System (DRPS) to generate the Experiment Data Records 
(EDR). • 

The IDR is recorded on a 9-track» 800 bpi, unlabeled magnetic tape. The 
recording method is NRZ-1 (non-retum-to-zero, change on ONE). A "I" 
bit is produced by reversal of flux polarity. The general format of 
the IDR file is shown in Figure A. 

IDR FILES 


An IDR file may be of any length, depending on the amount of data on the 
CODED DATA tape. It. shall consist of a File Label Record, Data Records, 
and an All Zeroes Record. An IDR file shall be an integral number of 
Data Records; if the last Data Record is not full, it is completed with 
filler blocks (defined below under Data Records). There shall be one 
IDR file on each output tape. The end of a tape is always indicated by 
a double end-of-file mark. 

IDR RECORDS 

1 

a. File Label Record . The File Label Record will always be the 
first record of an IDR file. It consists of 608 16-bit words 
(9728 bits) comprised of a 16-bit record sequence number (all 
zeroes), 8 zero bits, 5976 bits of data, and 3728 unused bits. 
The data portion will contain information specifying the IDR 
file type and contents. Table 1 describes the format of the 
IDR File Label Record. 

b. Data Records . The File Label Record is followed by N Data 
Records, where N i’s the’ appropriate number of records neces- 
sary to provide all data for a selected Probe, station, and 
pass. The format of the Data Record is shown in Table 2. 

Each Data Record consists of 608 16-bit words (9728 bits) com- 
prised of a 15-bit record sequence number (1 through 32767), 

7 additional 16-bit items of information, and 8 1200-bit 
blocks. The format of each 1200-bit block within the Data 
Record is shown in Table 3. Each block can be either a 
Telemetry Block or a Filler Block. (Filler Blocks are in- 
serted as the last blocks in an IDR File,' if necessary, to 
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complete an Integral number of Data Records.) A Telemetry 
Block consists of the Information shown in Table 2. A 
Filler Block consists of: 

Bits 1-48; same information described for 
the Telemetry Block in Tnble 2 

Bits 49-1200: repetitive 4210(g) pattern 

(i.e., 100010001000100010001000. . .) 

c. All Zeroes Record . The last record in each IDR file shall 
consist of; Bits 1-16: All I'a 

Bits 17-9728; All zeroes 


This record shall always precede the* double end-of-file 
mark which indlcatbat the end pf tape. 


FIGURE A - IDR FILE FORMAT 
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TABLE 1 - FILE UBEL RECORD FORMAT 


Bit No. 


1**16 

17-24 

25-48 

49-72 

73-96 

97-120 


Length 

(aits) 


16 

8 

24 

24 

24 

24 


121-152 

153-168 

169-192 

193-200 

201-216 

217-360 

361-384 

385-408 

409-432 

433-456 

457-9728 

Notes : 

(1) Right 

(2) Kanuali 


32 
16 
24 
8 

16 
144 
24 
24 
24 

24 

9272 

Ijustified 
Input. 


Content 


Description 


All zeroes 
All zeroes 
124 114 115 
.111 104 122 


( 8 ) 

( 8 ) 


031 


( 8 ) 


or 206 

^'(8) 

J^(8) 

30460 


( 8 ) 


( 8 ) 

30461 ( 3 ) 

All zeroes 

2 ASCII Characters 

■•FWD" or ''REV' 

®°(8) 

ni 

2 ASCII Characters'*^ 
All zeroes 
Binary Integer 
Binary Integer 
Binary Integer 


Binary Integer 


8 

9 

10 
n 


Record Sequence No. 

Spare 

3 ASCII Characters ("TLM") 

Tape ID; 3 ASCII Characters 
("IDR") 

DSS Number^^^’^^^ 

031 * DSS 14; 206 - OSS 43 


SP#1 

m wm 

SP#3 

LP4 


Probe ID 
(ASCII) 


( 2 ) 


( 2 ) 


Spare bits 
Tape Sequence No. 

Data Type; 3 ASCII Characters 
ASCII Character "0" 

Data Start Time, Year 
Bits not used 

Data Start Time, Day of Year 

Data Stop Time, Day of Year 

Data Start Time of Day (GMT) 
in 1/100 of second 

Data Stop Time of Day (GMT) 
in 1/100. of second 


( 2 ) 


All zeroes 


with leading character zeroes 
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TABLE 2 - DATA RECORD FORMAT 


Sit No. 


1-16 

17-32 
33-48 

49-64 
65-80 
81-128 
129-1328 
I 1329-2528 
2529-3728 
3729-4928 
4929-6128 
3129-7328 
7329-8528 
8529-9728 


Length 

(Bits) 


16 

16 

16 

16 

16 

48 

1200 

1200 

1200 

1200 

1200 

1200 

1200 

1200 


Content 


Binary Integer 


2300 

.0226 


( 8 ) 

( 8 ) 


0226 (gj 
0010 ( 8 ) 

All zeroes 
Data Block 
Data or Filler Block 


Note: 

(T)“lee Talile 3 for jlescription of each block. 


Description 


Record Sequence Nur.ber: 
Increments 1 for each Daia 
Record (from 1 to 32767) 

Record Size (in 8-bit bytes)' 

Size of each block (in 8-b1t 
bytes ) 

Same as bits 33-48 
Number of blocks in »*acord 
Spare bits 
1st Block^''^ 

2nd Block^^ '' 

3rd Block^^^ 

4th Block^"'^ 

5th Block 
5th Block 
7th Block 
8th Block 


(1) 

( 1 ) 

0 ) 

( 1 ) 







TABLE 3 » DATA BLOCK FORMAT 


Bit No. 

Length 
. (Bits) 

Content 

Descp1pt1on 

1-24 • 

24 

30473047 (Qj 

Sync Code 

25-32 

8 

031 ^gj oraoe^gj 

Identifies the OSS station whepe 
the data wepe received: 031 « 
OSSU; 206 - OSS 43(2) 

33-40 

8 

2=2(8) 

Destination Code 252 » ARC . 

41 -48 

8 

336(gj OP 276^gj 

Block Format Code: 336 * Data 
Block; 276 » Filler Block 
(Filler Blocks inserted as the 
last blocks of an IDR file, if 
necessary to complete an integral 
number of Data Records) 

49-51 ^ 

3 

®(8) 

Gross Data Description 

52-58^"*^ 

7 

0’3(8) 

User Dependent Type Code' 

59-65 ^ 

1 

ll4(gj OP 106(gj 

Data Dependent Type Code:® 114' ■ 
Large Probe; 106 » Small Probe 

66-72^''^ 

7 

or 011 ( 3 , 
or 012 ( 3 ) 
or 013(3, 

Spacecraft 



Small Probe #1 *010 
Small Probe 3*2-011 
Small Probe #3-012 
Large Probe - 013 

73-96^*’ ^ 

24 

Blnapy Integep 

Greenwich Mean Time (of day) in 
binary 1/100 seconds at which 
the last symbol of the block 
was received at the station 

97-98 


0 

Spare Bits 

99-108^'*^ 



Day of Year Code: 

(99-100)^' 
/ *! ' 

2 

Blnapy Integep 

Hundrisds digit 

(101-104)^' 

4 

Blnapy Integep 

Tens digit 

(105-108)^' 

4 

Blnapy Integep 

Units digit 

109-120^''^ 

12 

Blnapy Integep 

Block Serial Number which in- 
crements one count for aach 
data block in a file 


I 
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Table 3 - Da'ra fiiccv :ont'd.) 


Bit No. 

Length 

(Bits) 

Content 

121-128^’^ ■ 

8 

Binary Integer 

129-136^''^ 

8 

200 ( 3 , or 202 ( 3 ) 

137-160^''^ 

24 

All zeroes 

161-176^''^ 

16 

Two's Complement Binary 
No. 

177-192^^^ 

16 

Two's Complement Binary 
No. 

• 

193-704^^^ 

512 

Data for decoded frame or 
del eted frame 

■ 

705-720^^^ 

16 

Binary Integer 

721 -736^*'^ 

16 

1 

Binary Integer 


Description 


Millisecond Clock; a binary 
representation of the least 
significant decimal digit of 
the Greenwich Mean Time (de- 
rived from the same clock as 
that In bits 73-96) 

Telemetry Status: 200^gj « 
Pscoded Frame; 202 * Deleted 
Frame 

Spare Bits 

Symbol SNR^^^* bit 161 Is the 
sign bit (0»positive, 1 •nega- 
tive) and the binary point is 
located between bits 169 and 17C 

Receiver AGC^^^- bit 177 is ths 
sign bit (0*posit1ve, 1 “nega- 
tive) and the binary point is 
located between bits 185 and 18£ 


Decoded Frame: Bit 193 is the 

first bit of the frame and bit 
704 is the last bit of the 
Frame Sync Word (i.e., bits 681 
through 704 contain the Frame 
Sync Word; 11111000110001010100' 


001 


Deleted Frame: Bits 193 througt 

680 are all zeroes, bits 681 
through 704 contain the Frame 
Sync Word defined above. 


Number of symbol errors in the 
frame 


Number of computations perfcnnec 
in decoding the frame, in binary 
divided by 64, right justified 


) 
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Table 3 « Data Block Format (cont'd.) 


Bit Her. 


737-1200'’ 


Content 


464 All zeroes 


Description 


Unused bits 
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Enclosure 2 

P~V REVERSE PLAYBACK DECODING 
PRINTOUT FORMAT 


The computer printout shall contain the following information: 

1 . Title Block: 

a. Probe Identification 

b . DSS Number 

c. Data Start Time 

2. For Each Telemetry Frame Which is Printed : 

a. Frame number (i.e., block number) in decimal integer. 

b. Frame start time and stop time in hours, minutes, and 
seconds (seconds to 3 decimal places). 

c. Indication of whether it is a Deleted Frame. 

d. Indication of whether the frame was decoded by means of 
the "Quick-Look" decoding algorithm or the Fano algorithm; 
or whether it is a printout of the raw undecoded data 
symbols.’ 

e. Number of symbol errors corrected (where applicable). 

f. Average computations per b^’t, AC3 (where aoplicable). 

.The ACS is the number of computations performed, divided 

by the number of bits (512) in the frame. 

g. The 512 bits of decoded data or the raw data symbols 
(undecoded) in octal and binary form. 

3. Summary Block : 

a. Total number of frames processed. 

b. Total number of frames deleted (where applicable). 

c. Deletion rate (where applicable); i.e., number of frames 
deleted, divided by number of frames processed. 

d. Symbol error rate (where applicable); i.e., nu.mfaer of 
symbols corrected, divided by the total number of symbols 
processed. 
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